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CHAPTER  I 


INTRODUCTION 


1.1   Motivation  and  General  Purpose 

This  thesis  is  concerned  with  the  regulation  of  respir- 
atory muscles  in  animals.   The  living  organism  is  a  complex 
structure  with  many  degrees  of  freedom  and  whose  operation  is 
nonlinear  and  nons tationary .   The  deterministic  behavior,  if 
any,  is  usually  hidden  by  stochastic  variations.   The  ventila- 
tion system  of  any  living  organism  is  influenced  by  many  fac- 
tors.  An  excellent  conceptualization  was  suggested  by  Goodman 
(1964)  for  the  respiratory  system  of  man.   It  is  equally  ap- 
plicable to  any  living  organism  (Fig.  1.1). 
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Fig.  1.1.   Conceptual  diagram  showing  manipulation 
of  ventilation 


The  body  is  represented  as  an  active  system  containing  its 
own  store  of  internally  manipulated  energy.   Physical  processes 
are  coupled  with  the  central  nervous  system  and  the  environment. 
A  finite  number  of  state  variables,  e.g.,  blood  oxygen  tension, 
blood  CO^  content,  body  temperature,  etc.  indicated  by  X- ,  X„  , 
...,  X.  are  transformed  by  some  unidentified  dynamic  operators 
G..  ,  Gj  ,  ...,  G.  and  their  superimposed  effects  produce  the 
manipulation  of  ventilation.   Ideally,  knowledge  of  the  dynamic 
operators  G..  ,  G-,  ...,  G.  in  conjunction  with  the  state  of  the 
animal  would  enable  the  complete  prediction  of  the  organism's 
respiratory  behavior. 

The  purpose  of  this  thesis  evolved  from  research  currently 
underway  in  the  Department  of  Physiology  at  Kansas  State  Uni- 
versity.  Physiologists  are  investigating  chemoreceptors  asso- 
ciated with  the  C0_  sensing  mechanism  of  the  chicken.   These 
chemoreceptors  sense  the  blood  CO-  tension  and  stimulate  the 
neuro-muscular  system  to  adjust  the  respiratory  amplitude  to 
compensate  for  the  variation  in  blood  CO-  tension.   The  close 
relationship  between  a  dynamic  operator  of  Fig.  1.1  and  these 
chemoreceptors  is  obvious.   Experimental  data  consisted  of 
electrical  waveforms  recorded  from  within  the  chicken  respira- 
tory system  for  controlled  levels  of  blood  CO-  tension. 

Based  on  the  close  relationship  between  a  dynamic  operator 
and  the  C0„  chemoreceptors  it  was  asked  whether  a  model  ac- 
counting for  the  effect  of  the  CO.  chemoreceptors  on  respiratory 
electrical  activity  could  be  constructed.   Such  a, model  would 


lead  to  the  explicit  definition  of  at  least  one  dynamic  oper- 
ator in  Fig.  1.1  for  the  chicken.   Experience  and  knowledge 
obtained  from  the  chicken  (a  rather  primitive  organism)  would, 
it  is  hoped,  be  a  stepping  stone  to  the  understanding  of  the 
vital  process  of  respiration  and  body  functioning  in  higher 
animals,  including  man. 

This  thesis  presents  the  results  of  an  effort  to  illum- 
inate the  relationship  between  the  C0„  chemorecep tors  and 
respiratory  electrical  activity  through  an  analysis  of  the 
respiratory  electrical  activity  for  controlled  levels  of  CO- 
in  the  chicken.   In  addition,  a  mathematical  model  of  the 
respiratory  system  is  constructed  and  a  function  dependent 
on  only  the  system  parameters  and  state  of  the  chemoreceptors 
is  derived. 

It  is  hoped  that  the  results  of  this  work  will  be  carried 
on  until  all  the  dynamic  operators  associated  with  respiration 
are  discovered  and  the  conceptual  diagram  of  the  respiratory 
system  as  shown  in  Fig.  1.1  is  completely  described. 

1.2   Summary  of  Chapter  Development 

Chapter  II  begins  the  study  of  the  observed  C0_  chemore- 
ceptors effect  on  the  nerve  and  muscle  electrical  activity  with 
a  brief  introduction  to  the  physiological  processes  involved 
in  respiration.   Topics  include  the  generation  of  electrical 
behavior  and  the  observed  effect  of  CO-  on  the  respiratory 
system. 


A  model  of  the  respiratory  system  comprises  the  substance 
of  Chapter  III.   A  fairly  detailed  description  of  the  labora- 
tory experiment  is  presented  in  Chapter  IV.   A  block  diagram 
of  the  recording  system  is  also  included  in  Chapter  IV. 

Chapter  V  defines  the  mathematical  functions  that  are 
pertinent  in  studies  of  respiratory  electrical  behavior. 
The  cross  correlation  function  Is  shown  to  be  a  measure  of 
system  model  accuracy  in  this  chapter.   The  subject  of 
Chapter  VI  is  the  effect  of  digitizing  the  experiment  analog 
signals  and  the  effects  of  aliasing  due  to  a  sampling  rate 
below  the  Nyquist  frequency.   The  Fast  Fourier  Transform  is 
also  discussed. 

The  results  and  conclusions  based  on  computer  processing 
are  included  in  Chapter  VII.   In  addition,  at  the  back,  of  this 
thesis  are  included  two  appendices,  Appendix  A  is  a  description 
of  the  digital  processing  system  designed  for  the  data  pro- 
cessing required  by  this  thesis.   Appendix  B  is  a  mathematical 
analysis  that  is  used  to  derive  part  of  the  results  of  Sec. 
5.8. 


CHAPTER  II 


DESCRIPTION  OF  PHYSIOLOGICAL  PROCESSES 
INVOLVED  IN  RESPIRATION 


2.1   The  Nervous  System 

In  order  to  describe  the  physiological  processes  that  enter 
into  the  nerve  and  muscle  processes  leading  to  respiration,  a 
history  of  a  nerve  signal  will  be  traced  as  it  is  conducted 
along  a  fiber,  traverses  a  synapse  and  becomes  integrated  with 
other  electrical  activity,  and  finally  participating  in  the 
firing  of  a  muscle  fiber. 

2.1.1   Conduction  in  Fibers 

It  is  convenient  to  compare  neural  activity  to  an  electric 
circuit.   The  circuit  carries  messages  (nerve  signals)  through 
wires  (nerve  fibers)  to  a  load  (muscle  cells). 

A  major  difference  exists  between  the  nervous  system  and 
the  electrical  analogy.   A  nerve  fiber  has  the  ability  to 
propogate  a  signal  throughout  its  length  without  attenuation. 
This  is  possible  because  energy  is  available  to  the  nerve 
from  properties  associated  with  the  chemical  composition  of 
the  fluid  surrounding  and  internal  to  the  nerve  fiber.   Of 
course,  this  differs  from  a  wire  in  an  electrical  circuit  in 
which  a  transmitted  signal  is  attenuated  by  the  resistance  of 
the  wire. 


A  stimulus  produced  by  a  change  in  the  animal's  environ- 
ment will  cause  a  nerve  impulse  to  be  transmitted  along  the 
nerve.   If  the  stimulus  is  above  a  certain  threshold  level,  a 
pulse  will  be  conducted  along  a  fiber  (axon).   This  is  commonly 
referred  to  as  the  "all  or  none"  law. 


2.1.2   Electrical  Phenomena  Associated  with  Impulse 
Generator  and  Impulse  Conduction 


The  transmembrane  potential  of  a  nerve  cell  (neuron)  is 
due  to  the  separation  of  charge  across  the  membrane  surrounding 

the  neuroplasm  of  the  nerve.   The  membrane  of  a  resting  nerve 

+  + 

resists  passage  of  sodium  ions  (Na  ),  but  potassium  ions  (K  ) 

and  chloride  ions  (CI  )  move  more  freely  through  the  membrane. 

Outside  the  membrane  there  is  a  high  concentration  of  Na   and 

+  + 

CI   ions  and  a  low  concentration  of  K   ions,  since  the  K   ion 

can  move  out  of  the  cell  along  its  concentration  gradient  more 

easily  than  Na   ion  can  move  into  the  cell,  there  tends  to  be 

a  slight  excess  of  positively  charged  ions  in  the  outside  of 

the  membrane. 

When  a  resting  nerve  is  stimulated  above  the  threshold 

level,  a  wave  of  depolarization  travels  away  from  the  point  of 

stimulation.   This  wave  involves  an  increased  permeability  of 

Na   ions  through  the  membrane.   This  causes  the  inside  of  the 

cell  to  lose  its  normally  negative  charge  and  become  positive. 

Furthermore,  the  region  outside  the  membrane  having  lost  Na 

ions  now  is  at  a  lower  potential  than  the  area  still  surrounded 

by  Na   ions . 


Immediately  after  depolarization  occurs,  repolarization 
begins  as  the  result  of  an  increased  outward  flow  of  K   ions. 
This  reestablishes  the  separation  of  charge  and  hence  the  nor- 
mal resting  potential.   However,  the  interior  of  the  cell  now 
has  a  few  more  Na   ions  than  it  had  before  the  depolarization 
and  the  extra  cellular  fluid  now  has  a  few  more  K   ions.   Na 
ions  are  actively  transported  from  the  interior  of  the  nerve 
through  the  membrane  to  the  fluid  surrounding  the  nerve  and 
the  K   ions  are  actively  transported  to  the  interior  of  the 
cell.   The  original  concentration  gradient  is  then  restored 
across  the  membrane. 


2.1.3  Synaptic  Junctions 

The  neuron  is  composed  of  three  parts  1.   dendrites, 
2.   body,  and  3.  axon.   These  are  analogous  to  the  inputs,  net- 
work, and  outputs  of  an  electrical  circuit.   In  this  case  the 
network  is  a  simple  additive  type.   The  input  to  a  neuron  can 
consist  of  many  pulses  carried  by  axons  originating  from  other 
neurons.   The  output  of  the  neuron  is  the  integrated  effect  of 
the  input  impulses  and  is  directed  towards  other  nerve  cells 
or  toward  the  muscle  cells.   Such  a  junction  between  two  nerve 
cells  is  called  a  synaptic  junction. 

2.1.4  The  Nerve-Muscle  Junction 

The  end  of  a  nerve  leaving  the  spinal  cord  which  is  asso- 
ciated with  respiration  is  in  a  respiratory  muscle.   The 


electrical  signals  carried  by  the  axons  stop  here  and  do  not 
travel  beyond  the  nerve  ending.   Instead,  a  chemical  is  liber- 
ated at  the  interface  which  depolarizes  the  muscle  fibers.   This 
depolarization  in  turn  leads  to  contraction. 

2.1.5   Electrical  Action  of  Muscle  Cells 

When  a  muscle  fiber  is  stimulated  a  chemical  change  occurs 
over  its  entire  surface.  The  chemical  change  results  in  an  ex- 
ternal voltage  potential  in  the  region  surrounding  the  cell. 

2.2   Physical  Description  of  the  Respiratory  System 

The  respiration  system  in  the  chicken  differs  from  the 
respiratory  system  common  to  mammals.   In  the  mammal,  the  rate 
of  respiration  directly  affects  the  oxygen  (.0^)    and  carbon 
dioxide  (CCO  partial  pressures  in  the  blood.   A  similar  situa- 
tion exists  in  the  chicken.   However,  this  condition  can  be 
suppressed  surgically  without  drastically  altering  the  chicken's 
normal  body  functions. 

2.2.1   The  Respiratory  System 

Fig.  2.1  is  a  simplified  diagram  illustrating  the  basic 
action  and  parts  of  the  respiratory  system  of  the  chicken. 
The  air  sac  can  be  imagined  to  operate  as  a  fireplace  bellows. 
As  the  chest  muscles  squeeze  the  air  sac,  air  is  forced  out  of 
the  sac  and  through  the  air  ducts.  When  the  chest  muscles  ex- 
pand the  air  sac,  air  rushes  through  the  air  ducts  to  fill  the 
low  pressure  area  created  in  the  air  sac. 
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Fig.  2.1.   Simplified  respiratory  system 


As  muscles  cause  the  air  sac  to  expand,  air  is  drawn  into 
the  trachea.  In  the  animal,  the  number  of  parallel  air  ducts 
is  actually  much  greater  than  shown  in  the  figure.  Their  com- 
bined passage  area  is  much  greater  than  the  area  of  the  direct 
path  through  the  lung  to  the  air  sac  so  most  of  the  inrushing 
air  flows  through  the  passage  of  least  resistance,  namely,  the 
air  ducts. 

Gas  diffuses  into  the  diventricula  of  these  parallel  air 
ducts,  called  air  capillaries.   The  exchange  of  0„  and  CC^  oc- 
curs between  the  gas  and  the  blood  in  the  air  capillaries.   The 
efficiency  of  the  system  is  increased  by  the  unused  air  stored 
in  the  air  sac.   When  the  volume  of  the  body  decreases  on  ex- 
piration and  gas  is  forced  from  the  air  sac,  the  unused  gas 
located  therein  passes  once  wore  through  the  air  ducts  and  is 
available  for  exchange  a  second  time. 
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2.2.2   Modification  of  the  Respiratory  System  to  Enable 
Laboratory  Control  of  Respiratory  Movements 


A  mid-ventral  incision  from  the  chicken's  trachea  to  its 
cloaca  exposes  the  air  sacs.   The  air  sacs  are  then  punctured 
with  a  sharp  instrument.   Inserting  a  tube  into  the  trachea 
from  an  externally  controlled  gas  source  completes  the  opera- 
tion (Fig.  2.3).   Gas  is  then  passed  in  a  unidirectional 
manner  into  the  trachea,  through  the  lung,  and  out  to  the 
atmosphere  through  the  hole  in  the  air  sacs.   While  this  is 
done  the  animal  is  under  deep  anesthesia  and  has  no  sense  of 
pain.   If  exposed  tissues  are  kept  moist  with  physiological 
saline  the  bird  will  live  for  many  hours  in  this  state.   An 
excellent  description  of  this  surgical  technique  is  given  by 
DeWet  (1966,  pp.  250-263). 

As  the  C0„  chemoreceptor  in  the  animal's  body  senses  the 
C0„  concentration  of  the  gas  being  forced  through  its  air 
ducts  they  respond  by  stimulating  the  central  nervous  system  to 
manipulate  the  chest  muscles  in  a  manner  which  would  normally 
cause  the  C0„  level  to  return  a  normal  value.   With  the  air 
sacs  punctured  these  efforts  are  futile  because  the  bellows 
has  a  large  hole  in  it.   The  steady  stream  of  gas  which  is 
forced  through  the  air  ducts  is  unaffected  by  the  movement  of 
the  chest  muscles. 

To  increase  the  amplitude  of  respiration  (i.e.,  magnitude 
of  the  movement  of  the  chest  muscles)  it  is  only  necessary  to 
increase  the  percentage  of  C0„  in  the  gas  forced  through  the 
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Fig.  2.2.  Midventra.1  incision  in  Gallus  Domes ticus 
showing  position  of  recording  electrodes 
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lungs.   Decreasing  the  percentage  of  CO   in  the  ventilating 
gas  produces  the  opposite  reaction  and  respiratory  movement 
is  decreased.   It  is  even  possible  to  halt  respiratory  move- 
ment completely  by  eliminating  the  CO-  from  the  input  gas 
source . 
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CHAPTER  III 


A  MODEL  OF  THE  RESPIRATORY  SYSTEM 


3.1   A  Model  of  the  Respiratory  Nerve-Muscle  System 


From  the  previous  section  involving  the  physiology  of  the 
chicken,  a  model  representing  the  gross  behavior  of  the  respir- 
atory system  was  constructed.   It  includes  the  three  major  ele- 
ments of  the  system.   Namely,  the  CO _  sensing  dynamic  operator, 
the  central  nervous  system,  and  the  respiratory  muscles.   This 
model  is  shown  in  Fig.  3.1. 
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Fig.  3.1.   Model  of  the  chicken  respiratory  system 
including  feedback  and  instrumentation 
effects 


s'(t) 


3.2   Explanation  of  the  Model 


The  dynamic  operator  studied  is  represented  by  the  first 
box.   It  contains  the  unknown  dynamic  operator  that  senses  the 
inspired  CO.  level  and  triggers  the  nervous  system  to  respond 


u 


appropriately.   The  central  nervous  system  and  its  myriad  of 
nerve  bundles  is  represented  by  the  second  box.   It  is  assumed 
to  contain  only  those  nerves  associated  with  the  respiratory 
system.   The  electrical  signals  carried  by  all  the  nerves  are 
called  collectively  r(t). 

The  bird  responds  to  the  electrical  signals  of  the  nervous 
system  by  increasing  or  decreasing  the  contraction  of  its 
respiratory  muscles.   The  muscles  involved  in  this  movement 
are  represented  by  the  third  box.   The  muscles  react  to  the 
stimulating  electrical  signals  (Refer  to  Sec.  2.1.5)  which  are 
represented  by  s(t)  in  the  model. 

Ordinarily,  the  movement  of  the  chest  muscles  causes  an 
increase  or  decrease  of  C0„  in  the  animal's  blood.   This  causes 
the  C0_  sensing  dynamic  operator  to  change  its  state.   This 
action  is  represented  by  the  feedback  loop  from  the  muscle 
output  to  the  CO-  sensing  chemoreceptors  in  Fig.  3.1. 

F..  and  F„  represent  the  bandwidth  limitations  imposed  on, 
the  electrical  activity  by  the  recording  instrumentation.   The 
signals  perturbed  by  instrumentation  are  called  r'(t)  and  s  '  (t) 
for  the  perturbed  nerve  and  muscle  signals,  respectively. 


3.3   Limitation  of  the  Model  Associated  with  Constraints 
Imposed  by  Instrumentation 


It  is  apparent  that  the  compactness  and  delicate  nature  of 
the  electrical  probes  used  to  record  the  animal's  electrical 
activity  limits  the  application  of  the  model  to  a  discrete 
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level  of  the  overall  respiratory  system,  mainly  due  to  the 
fact  that  only  the  electrical  activity  of  one  nerve  bundle 
can  be  recorded.   Likewise,  only  the  muscular  activity  from 
a  small  area  near  the  recording  electrode  was  recorded.   It 
is  hoped  that  the  results  of  this  sampling  are  applicable  to 
generalizations  about  the  entire  system.   However,  this  can- 
not be  verified  at  this  time. 
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CHAPTER  IV 


DESCRIPTION  OF  THE  EXPERIMENT  PERFORMED 
ON  THE  CALLUS  DOMESTICUS 


4.1  Introduction 

Observation  of  the  respiratory  movements  in  the  Gallus 
Domes  tlcus  (common  chicken)  shows  that  the  respiratory  muscles 
undergo  a  smooth,  graded  contraction.   This  is  brought  about 
by  a  host  of  complex  properties,  for  example,  muscle  structure, 
membrane  properties,  patterns  of  innervation,  and  excitation- 
contraction  coupling  mechanisms.   This  section  illustrates  the 
techniques  used  to  record  electrical  activity  associated  with 
the  relationship  of  the  C0„  dynamic  operator's  state  upon  the 
electrical  waveforms  in  the  respiratory  nerves  and  muscles 
under  its  control. 

4.2  Method  of  Varying  Nerve  and  Muscle  Electrical 
Activity  with  Selected  Levels  of  C0» 

The  phenomenon  studied  was  the  influence  of  the  dynamic 
operator  as  it  manifested  itself  in  the  electrical  waveforms 
of  the  respiratory  nerve  and  muscle  system  for  various  levels 
of  C02.   Thus  it  was  first  necessary  to  prepare  the  animal  for 
external  control  of  the  dynamic  operator  and  respiratory  rate 
as  described  in  Sec.  2.2.2.   The  gas  mixture  forced  through  the 
trachea  consisted  of  air  supplied  by  a  pressure  vacuum  air  pump 
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and  C0„  supplied  from  a  tank.   These  gases  were  metered  with  a 
floating  ball  meter  for  air  and  a  differential  pressure  meter 
for  CO,,.   Since  air  contains  very  little  CO-,  the  C0„  in  the 
animals  lungs  was  controlled  by  regulating  the  C0„  from  the 
tank. 

The  percentage  of  C0_  varied  the  state  of  the  respiratory 
muscles.   When  the  percentage  was  low  (below  3-6%)  the  animal 
ceased  phasic  contractions  and  exhibited  no  respiratory  move- 
ment (state  of  apnea).   When  the  percentage  was  elevated  the 
animal  began  phasic  respiratory  movements  which  became  pro- 
gressively greater  as  the  C0_  level  was  increased.   The  levels 
of  C02  used  were  0%,  2%,  4%,  6%,  8%,  10%,  12%,  and  14%. 

4.3   Electrical  Monitoring  and  Recording  Apparatus 

The  recording  equipment  was  substantially  the  same  as  used 
by  De  Wet  (De  Wet,  1966).   The  only  modification  was  the  addi- 
tion of  a  Precision  Instruments  FM  one  channel  tape  recorder 
(Fig.  4.1).   This  recorder  provided  a  time  expansion-contraction 
capability.   It  was  possible  to  record  and  playback  at  37.5 
ips  and  .375  ips  respectively,  without  loss  of  frequency  re- 
sponse . 

Electrical  recordings  of  efferent  neural  activity  to  a 
respiratory  muscle  were  taken  from  a  tiny  isolated  nerve 
branch  of  the  6th  thoracic  nerve.   The  nerve  was  cleaned  of 
connective  tissue  and  draped  over  a  small  pair  of  silver  elec- 
trodes.  The  electrodes  were  attached  to  a  cathode  follower 
impedance  matching  amplifier  and  then  into  a  high  gain  ac 
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amplifier.   The  amplifier  output  was  parallel  to  one  channel  of 
a  multichannel  oscilloscope  and  to  the  normally  closed  position 
of  a  coax  switch.   The  pole  of  the  coax  switch  was  connected 
to  the  FM  recorder  input  (See  Fig.  4.1). 

Electrical  activity  from  the  respiratory  muscle  (M. 
Transversus  Abdominis)  was  recorded  with  two  small,  flexible, 
teflon-coated,  stainless  steel  electrodes.   These  electrodes 
were  sharpened  so  they  easily  penetrated  the  muscle.   The 
electrodes  were  attached  to  another  impedance  matching  cathode 
follower  and  to  another  high  gain  ac  amplifier.   The  amplifier 
output  was  parallel  to  the  oscilloscope  and  to  the  normally 
open  position  on  the  coax  switch. 

Energizing  the  coax  switch,  switched  the  recorder  input 
from  the  nerve  electrodes  to  the  muscle  electrodes.   The  entire 
system  was  calibrated  so  that  80  pvolts  at  the  electrodes 
corresponded  to  1.4  volts  output  on  the  FM  tape  recorder.   The 
entire  recording  system  was  located  outside  the  wire  screened 
operating  room  to  minimize  electrical  noise  from  60  Hz  sources. 

4.4   Recording  Procedure 

It  was  found  that  high  levels  of  C0„  stimulated  the  animal 
to  the  extent  that  his  movements  were  great  enough  to  tear  out 
the  carefully  implanted  electrodes.   Thus  it  was  decided  to 
begin  recording  nerve  and  muscle  activity  for  high  C09  levels 
and  then  decrease  to  lower  levels.   This  allowed  the  critical 
stimulation  point  to  occur  when  the  animal  was  under  maximum 
anesthesia. 
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Fig.  4.1.  FM  recorder  and  coax  switch  used  to  record 
nerve  and  muscle  electrical  activity 
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Ten  seconds  of  nerve  electrical  activity  were  recorded 
immediately  followed  by  ten  seconds  of  muscle  electrical 
activity  for  C02  levels  of  14%,  12%,  10%,  8%,  6%,  4%,  2%,  and 
0%  in  that  order.   After  the  C0„  level  was  changed  each  time, 
a  delay  of  three  minutes  before  recording  was  maintained  to 
let  the  animal  adjust  completely  to  the  new  C02  level. 

4.5   Brief  Description  of  Data  Processing  System 

Appendix  a  contains  a  detailed  discussion  and  explanation 
of  the  analog-to-digital  data  processing  system  designed  and 
developed  for  use  in  biomedical  research  at  Kansas  State  Uni- 
versity.  The  following  sections  serve  to  illuminate  the 
process  applied  to  the  experimental  data,  recorded  as  pre- 
viously described  in  order  for  it  to  be  analyzed  by  a  digital 
computer. 

Basically,  the  processing  system  samples  and  quantizes 
analog  data  into  a  digital  format  suitable  for  computer  in- 
put.  The  maximum  sample  and  quanti tization  rate  is  50  samples 
per  second.   This  rate  is  usually  sufficient  for  most  bio- 
medical data  (e.g.,  electrocardiograms,  fetal  heartbeats, 
finger  twitch,  etc.).   Higher  sample  and  quantizing  rates  can 
be  achieved  by  time  expansion  techniques  (Refer  to  Sec.  6.1). 

The  data  processing  system  consisted  of  an  analog-to- 
digital  convertor  (ADC)  and  a  digital  stepping  recorder  (DSR) 
both  under  the  control  of  a  small  general  purpose  computer. 
The  ADC  was  manufactured  by  Texas  Instruments  and'is  discussed 
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in  Sec.  A. 2.       Digi-Data  Corporation  produced  the  DSR  which  is 
discussed  in  Sec.  A. 3.   The  general  purpose  computer  was  built 
by  Remington-Rand.   The  computer  and  the  machine  language 
programs  used  in  the  control  process  can  be  found  in  Sec.  A. 5. 

4.5.1   Creation  of  a  Digital  Magnetic  Tape  Containing 
Biomedical  Data 

The  general  purpose  computer  (Athena)  controlled  the 
sampling  and  recording  equipment  only.   It  performed  no  data 
analysis.   Analysis  was  done  later  with  the  IBM  360/50  com- 
puter using  the  magnetic  tape  generated  by  the  Athena  con- 
trolled processing  system.   Operation  of  the  sampling  and 
quantitization  system  began  with  Athena  addressing  a  selected 
analog  input  channel  on  the  ADC.   Quantitization  of  a  40  usee, 
portion  of  the  analog  signal  was  done  automatically  by  the 
ADC  1.5  usee,  after  receiving  the  channel  address  command  from 
the  computer. 

The  14  bit  ADC  output,  in  offset  binary  code,  was  then 
stored  in  the  256  twenty-four  bit  word  core  of  the  computer. 
The  computer  then  placed  the  digitized  sample  and  channel 
address  at  the  DSR  input  where  it  was  recorded  on  magnetic 
tape  in  a  specified  code  (See  Fig.  A. 8).   This  sequence  oc- 
curred every  5  ysec.  ,  yielding  a  real-time  sampling  rate  of 
50  samples  per  second. 
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4.5.2   Decoding  of  the  Coded  Biomedical  Data  by  the 
IBM  360/50  Computer 


The  magnetic  tape  produced  by  the  data  processing  system 
containing  experimental  data  was  not  in  a  standard  IBM  code 
(e.g.,  BCD,  EBCIDEC,  etc.).   Therefore,  further  processing  was 
required  on  the  data  before  a  useful  analysis  could  be  made. 

The  magnetic  tape  units  associated  with  the  IBM  360/50 
have  a  capability  of  reading  binary  data  into  core  without 
translation.   Fig.  A. 9  (Appendix  A)  illustrates  the  way  eight 
seven  track  tape  characters  were  located  in  core  bytes  after 
such  a  read  occurred  on  binary  data  without  translation.   Fig. 
A. 11  is  a  flow  chart  of  an  assembler  language  program  that 
converted  the  binary  data  in  core  back  to  a  form  differing  from 
the  measured  voltage  by  only  a  scale  factor.   The  assembler 
program  writes  the  converted  data  on  a  nine  track  tape  unit. 
The  nine-track  tape  was  then  called  by  a  standard  Fortran  IV 
program  which  multiplied  by  the  scale  factor  and  located  the 
data  once  more  in  core  but  with  the  proper  format  for  computer 
analysis . 
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CHAPTER  V 


DEFINING  THE  MATHEMATICAL  FUNCTIONS  PERTINENT  TO 
THE  ANALYSIS  OF  THE  RECORDED  ELECTRICAL  ACTIVITY 


5.1   Introduction 

The  mathematical  functions  that  were  thought  to  be  signi- 
ficant in  the  study  of  nerve  and  muscle  electrical  activity  are 
defined  in  this  chpater.   Although  not  all  were  used  in  this 
thesis,  it  is  felt  that  defining  them  here  will  serve  to  benefit 
future  work  in  this  area.   An  analysis  of  a  mathematical  model 
of  the  respiratory  system  is  presented  in  Sec.  5.8.1.   The 
mathematical  function  used  mainly  in  the  subsequent  analysis 
of  the  electrical  activity  was  the  Fourier  transform  and  the 
associated  power  amplitude  spectrum.   The  reason  for  not  using 
the  remaining  functions  was  two  fold.   One,  computer  computation 
time  exceeded  reasonable  limits.   Two,  the  significance  of  the 
crosscorrelation  function  for  two  waveforms  not  recorded  simul- 
taneously yields  little  information  of  value  if  their  similarity 
at  the  sair.e  time  is  to  be  studied.   The  mathematical  functions 
of  interest  are: 

1.  the  Fourier  Series  Amplitude-Frequency  Plot 

2.  the  Fourier  Series  Phase-Frequency  Plot 

3.  the  Fourier  Integral 

4.  the  Crosscorrelation  Functi  i 

5.  the  Fourier  Frequency  Plots  of  the  Crosscorrelation 
Function . 
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5.2   The  Fourier  Series 

If  a  function  f(t)  is  a  periodic  function  of  time  with 
period  T,  it  may  be  expressed  as  a  Fourier  Series  as  in  Schwartz 
(1959,  pp.  21-23)  by, 


f (t)  -  —  +  -   I  (&      cosw  t  +  b   sinw  t) 
1    T    .   n      n      n      n 
n=l 


(5.1) 


2rrn 
where  ,      oi       •*    ~~^~~ 


The  coefficients  a   and  b   are  found  through  the  equations, 

n       n  °  n 


T/2 
a   =  /     f(t)  cosu  t  dt    n=0 ,  1.  ... 
n    -T/2 


(5.2a) 


T/2 
b   ■  /     f(t)  sinw  t  dt    n=l,  2,  ... 
n    -T/2 


(5.2b) 


5.3   The  Amplitude-Frequency  Plot 


The  amplitude-frequency  plot  is  proportional  to  a  plot  of 


/  2   2   versus  cj  .   Often  this  olot  is  referred  to  as  the  ampli- 

a  +b  n  '  r 

n   n 

tude  spectrum  of  the  function  f(t). 


5.4   The  Phase-Frequency  Characteristic  Plot 


The  phase-frequency  characteristic  plot  is  a  plot  of 


tan  (-b  /a  )  versus  to 
n   n  n 
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5.5   Complex  Exponential  Fourier  Series  Form 

This  is  an  alternative  form  of  the  Fourier  series  written 


as  , 


f  (t)    =   -4  +   V  ■    S    C       cos(w    t+0    ) 
•lvw  TT  n  nn 

n="l 


(5.3) 


where,         Cn    =    ^2+b2 

n      n 


0      =    tan^C-b    /a    ) 
n  n      n 


u 


2irn 


n 


5.6      The    Fourier    Transform 


If    a    function    is    not    periodic    the    Fourier    series    definition 
of    Eq.    5.1    is    not    applicable.       Examples    of   nonperiodic    functions 
are    transient    responses    and    unit    impulse    responses.       A    Fourier 
transform    can    be    defined    for    this    case    by    the    equation, 


F(ju)    =    /    f(t)e"jU)tdt 


(5.4a) 


and  the  corresponding  time  function  can  be  recovered  from  F(jo>) 
with  the  equation, 


f  (t)  =  -~  f    F(jw)eja)tdu> 


(5.4b) 
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provided,      /  |f(t)|dt   is  finite. 


5.7   The  Crosscorrelation  Function 

If  f,  (t)  and  f2(t)  are  different  periodic  functions  of  the 
same  fundamental  frequency,  the  crosscorrelation  function  of 
these  functions  is  defined  in  Lee  (1960,  pp.  17-18)  as, 

,   T/2 
Q19(T)  -£/    f,(t)  f  (t+T)dt  (5.5a) 

11  £    -T/2  X 

where,  T  is  the  period  of  f,(t)  and  f2(t) 
T  is  the  time  displacement. 
For  aperiodic  functions,  the  crosscorrelation  function  is 
defined  as , 

CO 

*,9(T)  =   /  f,(t)  f2(t+T)dt  (5.5b) 

where  f.,(t)  and  or  f2(t)  are  aperiodic. 

5.8   Fourier  Series  of  the  Crosscorrelation  Function 
Writing  the  Fourier  series  of  f^t)  and  f2(t)  of  period  T 


as  , 


(c)  =  _i£-  +  1   E  C,   cos(w  t+9,  )  (5.6) 

v  '  T     T    ..  In       n    In 

n  =  l 


f  (t)  -  -~  +  |  2  C2a  cos(unt+62a)  (5.7) 

n=l 
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2irn 
and      to      =    — — 
n  T 


enables    the    Fourier    series    of    the    crosscorrelation    function    to 
take    the    form, 

a   a  °° 

*,,(T)  =   1092°  +  ~   E  C,  C,   cos(o)  T+60  -9.  )        (5.8) 
12         „2      _2      In  2n       n    2n   In 
1       i   n=  1 


where,    C.   =  /  2  .,  2 
'     in     a .  +b  . 
jn   jn 


I.   =  tan_1(-b.  /a.  )   .  (5.9) 

]"  jn   jn 


5.8.1   The  Amplitude-Frequency  Plot  of  the  Crosscorrelation 
Function 


The  amplitude-frequency  plot  of  $.2(t)  is  a  plot  of  C    C, 

versus  to  .   In  terms  of  the  original  time  functions  this  is  a 

n  ° 


plot  of  /  2  , ,  2  •  /l     ,'  2   versus  oj  .   As  a  consequence,  if  the 
a- +b,     a0  +b0  n 

la  In     2n   2n 

average  term  or  any  harmonic  is  absent  in  either  f.(t)  or  fjCt), 
the  constant  term  or  the  corresponding  harmonic  term  will  be 
absent  in  the  Fourier  series  of  the  crosscorrelation  function. 


5.8.2   The  Phase-Frequency  Plot  of  the  Crosscorrelation 
Function 


The  phase-frequency  plot  of  the  crosscorrelation  function 

is  a  clot  of  9_  -6,   versus  u  .   In  terms  of  the  original 
2n   In  n  a 

functions  this  is  a  plot  of  the  difference  in  phase  between 

f.  (t)  and  f0(t)  versus  to  . 
I  l  n 
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5.8.3   A  Mathematical  Analysis  of  the  Respiratory  System 
Under  Open  Loop  Conditions 


The  purpose  of  the  following  analysis  is  to  demonstrate  the 
applicability  of  the  previously  defined  mathematical  functions 
to  the  respiratory  system  model  of  Sec.  3.1  when  the  feedback 
loop  is  broken  as  described  in  Sec.  2.2.2.   The  analysis  pre- 
sented is  a  calculation  of  the  crosscorrelation  function  between 
the  nerve  and  muscle  electrical  waveforms.   The  crosscorrelation 
function  was  chosen  because  both  the  nerve  and  muscle  signals 
are  available  simultaneously  to  the  researcher  so  that  the 
nerve-muscle  crosscorrelation  function  computed  with  experimental 
data  can  be  used  as  a  check  against  the  same  crosscorrelation 
function  derived  from  a  mathematical  model  of  the  respiratory 
system. 

The  following  analysis  yields  an  equation  relating  the 
block  diagram  components  h   (t),  h2  (t),  and  fpp(t)  to  the 
nerve-muscle  crosscorrelation  function.   Interpretation  of  this 
equation  reveals  that  the  crosscorrelation  function  is  suitable 
for  use  as  a  gauge  to  the  accuracy  of  respiratory  system  models. 

In  Fig.  5.1  the  signal  received  by  the  central  nervous 
system  from  the  C09  dynamic  operator  is  labeled  fpp(t).   The 
affect  of  the  central  nervous  system  on  fpp(t)  is  represented  by 
h.(t),  where  h^t)  is  the  linear  impulse  response  of  the  central 
nervous  system.   Thenerve  electrical  activity  s(t)  is  shaped 
by  h2(t)  to  produce  the  muscle  electrical  activity  r(t).   It  is 
assumed  that  h„(t)  is  the  linear  impulse  response  of  the  muscle 
sys  tem. 
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Fig.  5.1.   Mathematical  model  of  the  respiratory  system 
under  open  loop  conditions 


5.8.4   Derivation  of  the  Nerve-Muscle  Cross  correlation 
Function 


Due  to  the  phasic  nature  of  the  respiratory  electrical 
activity  the  cross  correlation  defining  equation  for  periodic 
functions  is  used  in  this  analysis.   The  crosscorrelation  be- 
tween r(t)  and  s(t),  both  of  period  T  is  written  as 


^c<T>  "  9T  f     r<C>   S(t+T)dt 

rs       2T  _T 


(5.10) 


where  $   (x)  is  the  nerve-muscle  crosscorrelation  function, 
rs 

Because  the  nerve  and  muscle  impulse  response  functions 
h1(t)  and  h„(t)  were  assumed  to  be  linear,  r(t)  and  s(t)  can  be 
written  in  terms  of  the  following  convolution  integrals 


r(t)  =  /  h2(v)  s(t-v)dv   , 

—  CO 


(5.11) 


and , 


s(t)  =  /  h1  (a)fpp  (t-a)dcr   , 

—  CO 


(5.12) 


where  v  and  G  are  dummy  variables.   Direct  substitution  of  Eqs. 
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5.11  and  5.12  into  Eq .  5.10  yields  for  the  crosscorrelation 
function , 


$   (t)-~  /  /  h  (v)s(t-v)dv/  h  (a)fPP(t+T-a)dodt   .  (5.13) 

—  "£   —CO  —  00 


Defining , 


*sf (T+v-a)=y-/  s(t-v)fpp(t+x-a)dt   , 


(5.14) 


Eq .  5.13  becomes, 


*   (x)  =  /   /  h0  (v)h..  (a)<2>  ,(t+v-  a)  dvda 
rs  I  1     si 

—  00   —00 


(5.15) 


Eq.  5.16  can  be  manipulated  so  that, 


(tW  h  (v)  [/  h..(a)4>   (x+v-a)da]dv   , 


rs  '       2    L    1     s  f 

—  00  —00 


(5.16) 


where  the  quantity  in  brackets  is  recognized  as  the  convolution 
of  the  impulse  response  h..  (t)  with  the  crosscorrelation  function 
of  the  CC_  operator  waveform  fpp(t)  and  the  nerve  waveform  s (t) . 
Using  the  symbol  *  to  indicate  convolution,  Eq.  5.16  can  be 
written  as , 


0   (t)  =  /  h0(v)h..  (a)*$  r(x+v-0)dv   . 
rs  I  1      si 


(5.17) 
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The  variables  associated  with  the  convolution  result  in  Eq.  5.17 
are  T  and  v.   With  this  fact  Eq.  5.17  is  recognized  as  the 
aperiodic  form  of  the  cross  correlation  function  of  the  impulse 
response  h0  (t)  and  the  convolution  result  indicated  in  Eq.  5.17. 
Therefore  Eq.  5.17  can  be  expressed  as, 


>rs(x)  -  h2(t)  e  rh]L(t)*4>sf(T)]    , 


(5.18) 


where  the  symbol  G  means  aperiodic  correlation  and  *  indicates 
convolution 

This  equation  can  be  expanded  further  by  writing  §  f(t) 
in  terms  of  the  block  diagram  components.   Pursuing  this  course, 

<*>   (x)  by  definition  is, 

rs 


$fs(x)  -  2T  /  s(t)fPP(t+T)dt 


(5.19) 


Substitution  of  Eq.  5.12  into  Eq.  5.19  for  s(t)  results  in 
the  equation , 


X     00 


°sf(T)  "  2T  /  f   h1(a)fPP(t-a)d<*fpp(t+T)dt   ,         (5.20) 


—  X  —oo 


which  can  be  reduced  easily  to, 


$  f(0  -  /  h  (a)*,,(T+a)da   , 

si  1   •  f f 


(5.21) 
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where , 


•ff(T) 


— -  /  fpp(t)fpp(t+T)dt   . 
-T 


In  terms  of  the  previous  notation  for  aperiodic  correlation, 
Eq.  5.21  becomes, 


»   (t)  =  h.  (t)  0  4>   (t)   . 

st        1        ft 


(5.22) 


Upon  substitution  of  this  result  into  Eq.  5.18,  the  cross- 
correlation  function  of  the  nerve  and  muscle  electrical  activity 
is  found  to  be, 


»rs(T)  =  h2(t)  0  {h1(t)*[h1(t)  0  *ff(T)]}   . 


(5.23) 


5.8.5   Comments  on  the  Significance  of  the  Nerve  Muscle 

Crosscorrelation  Function  in  Light  of  the  Mathematical 
Analysis 


It  can  be  concluded  from  Eq.  5.23  that  the  crosscorrelation 
function  of  the  nerve  and  muscle  waveforms  is  solely  dependent 
on  the  mathematical  model  system  parameters  and  the  autocorrelation 
function  of  the  dynamic  operator  waveforms. 

The  operations  called  for  by  Eq.  5.23  are  difficult  to  eval- 
uate for  any  choice  of  h,(t),  h„(t),  or  §ff(x)  even  with  a  large 
digital  computer.   With  the  advent  of  the  Fast  Fourier  transform 
(to  be  discussed)  it  is  a  relatively  easy  matter  to  compute  the 
Fourier  transform  of  Eq.  5.23  on  a  digital  computer.   Using  the 
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results  of  Appendix  B,  the  Fourier  transform  of  $   (t)  is 

rs 


*   (jw)  -  H_(-jw)|H.  (jw)|2  <J>_(jw)   , 


rs 


ff 


(5.24) 


where  H0  (-jio)  is  the  complex  conjugate  of  the  Fourier  transform 
of  h0(t),  |H„(ju))|  is  the  magnitude  of  the  Fourier  transform  of 
h.(t)  and  <I>  .(jo))  is  the  power  amplitude  spectrum  of  the  C0„ 
dynamic  operator  waveform. 

Thus  it  is  concluded  that  modeling  in  the  frequency  domain 
of  the  respiratory  system  is  the  best  method  at  present.   Eq. 
5.24  can  be  utilized  as  a  guide  to  the  accuracy  of  the  model 
since  the  cross  correlation  function  and  its  Fourier  transform 
can  be  computed  for  real  data  and  then  compared  to  the  results 
of  Eq.  5.24. 

Results  of  this  mathematical  analysis  indicated  that  the 
study  of  the  respiratory  system  begin  in  the  direction  of  un- 
covering its  characteristics  in  the  frequency  domain.   This  is 
the  course  of  investigation  that  is  followed  throughout  the 
remainder  of  this  thesis. 
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CHAPTER  VI 

METHOD  OF  COMPUTER  ANALYSIS 

6.1   Introduction 

Two  animals  were  used  in  the  experiment  from  which  nerve 
and  muscle  activity  was  recorded  for  each  level  of  C0»  as 
stated  in  Sec.  4.4.   Recording  was  done  with  a  Precision 
Instruments  FM  recorder  at  a  speed  of  37.5  ips  with  a  flat 
bandpass  from  OHz  to  10  KHz.   Selected  protions  of  each  sig- 
nal at  each  C02  level  were  converted  into  digital  form 
suitable  for  IBM  360/50  input  by  the  data  conversion  system 
described  in  Appendix  A.   The  data  conversion  system  sampled 
the  selected  portions  of  the  electrical  signals  50  times/sec. 
To  increase  the  effective  sampling  rate  of  the  waveforms  the 
FM  recorder  was  slowed  to  a  speed  of  .375  ips  and  its  band- 
width reduced  to  OHz  to  100Hz  before  processing  by  the  data 
conversion  system.   This  yielded  an  effective  sampling  rate  on 
the  original  signals  of  5000  samples/sec.   The  Nyquist  sampling 
theorem  states  that  bandlimited  signals  of  finite  energy  with 
bandwidth  f   may  be  completely  described  from  a  knowledge  of 
its  samples  taken  at  a  rate  of  2f   per  second  (Landau,  1967). 
The  highest  possible  frequency  in  the  waveforms  was  10,000  Hz 
due  to  the  recorder  bandpass  filter.   Thus  a  Nyquist  sampling 
rate  of  20,000  samples  would  be  required  for  complete  recovery. 
Since  the  sampling  was  done  at  a  rate  one-fourth  as  fast  as 
the  Nyquist  frequency,  there  exists  the  possibility  of  aliasing 
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from  high  frequency  components.   This  is  discussed  in  the 
next  section. 

6.2   The  Effect  of  Aliasing  Due  to  a  Sampling  Rate  Less 
That  the  Nyquist  Frequency 

A  discussion  of  aliasing  is  given  in  Blackman  and  Tukey 
(195S)  for  a  stationary  random  process  which  has  been  sampled 
at  uniformly  spaced  intervals  of  time  t  viz, 

t  =  0,  At,  2At,  ... 

Assuming  that  the  signals  in  this  thesis  are  at  least  stationary 
over  the  period  of  observation,  a  sampled  waveform  can  be  de- 
fined as 

f  (t)  =  V(t;At)  f(t)   ,  (6.1) 

s 

where  f  (t)  is  the  sampled  signal  obtained  from  the  analog 

signal  f(t)  by  equispaced  samples  At  sees,  apart  and  V(t;At) 

represents  an  infinite  dirac  comb  (Blackman  and  Tukey,  1958, 

p.  6  8)  where, 

n=°° 
V(t;At)  =  At  I    6(t-nAt) 

The  Fourier  transform  of  the  analog  signal  is 

,   T/2 
F(f)  =  i  /    f(t)  exp(-jut)dt   ,  (6.2) 

1  -T/2 

where    F(f)    is    the    real    Fourier    transform    and    u    =    2irf. 
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The  power  amplitude  spectrum  corresponding  to  F(f)  is 
simply , 

P(f)  =  |F(f)!2   ,  (6.3) 

where  P(f)  is  the  power  spectrum  of  the  original  time  function 
f(t). 

The  Fourier  transform  of  the  sampled  signal  is  given  by 

,   T/2 

F  (f)  -  ~  f         f  (t)  exp(-ju>t)dt   ,  (6.4) 

s       i  _T/2  s 

where  F(f)  represents  the  Fourier  transform  of  the  discrete 
sampled  function  f(t).   Substituting  6.1  into  6.2  yields 

,   T/2 
F  (f)  =  ~   /    V(t;At)f(t)  ex0(-jut)dt   .  (6.5) 

S       l    -T/2 

It  is  well  known  that  the  Fourier  transform  of  a  product  can 
be  obtained  by  convoluting  the  individual  Fourier  transforms 
in  the  frequency  domain.   Applying  this  fact  to  Eq.  6.5  yields 
the  Fourier  transform 

Fs(f)  =  A(f~)  *  P(f)   ,  (6.6) 

where  A(f;— )  is  an  infinite  dirac  comb  in  the  frequency  domain 
given  by, 

A(f;^)  =  "^(f-^)    , 
and  *  means  convolution. 


Eq.     6.6    can    be    expanded    to, 
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Fs(f)  -    z    nt-~) 

n  =  -°° 


E  F(f-^);n^o 


7\=>-< 


(6.7) 


Note  that  F  (f)  differs  from  F(f)  by  the  terms, 


Z       F(f-~)  ;  n^O 
n=-°° 


which  are  identical  in  shape  to  the  actual  Fourier  transform, 

and  centered  on  n/At.   The  effect  of  these  additional  terms 

on  F  (f)  is  called  aliasing  since  without  them  F  (f)  would  be 
s  s 

identically  F(f)  and  hence  P  (f)  would  equal  P(f). 

Recognizing  that  1/At  is  equal  to  the  sampling  frequency 

called  f  ,  then 
s  ' 

fs  -  5000  Hz 

for  the  data  processing  system. 
Eq .  6.7  becomes , 


F  (f)  =  F(f)  +  Z       F(f-nf  );  n^O 
s  s 

n=_co 

for  the  data  processing  system.   A  little  thought  reveals  that 

F  (f)  is  identical  to  F(f)  between  0  and  f  /2  Hz  if  the  condi- 
s  s 

tion  that  no  frequency  components  exist  in  F(f)  above  f  /2  is 
satisfied.   Therefore,  the  data  processing  system  can  be  ex- 
pected to  furnish  a  valid  representation  of  the  actual  power 
amplitude  plot  up  to  f  /2  or  2500  Hz  only  if  there  are  no 
frequency  components  in  the  experimental  waveforms  above  2500  Hz 
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The  components  in  the  recorded  waveforms  above  2500  Hz 
are  small  compared  with  the  lower  frequency  components.   A 
quantitative  demonstration  of  this  was  done  by  filtering  the 
experimental  waveforms  with  high  frequency  cutoffs  of  5000  Hz 
and  2500  Hz  and  observing  the  affects  on  the  waveform.   Only 
a  slight  affect  on  the  waveform  was  observed  for  the  5000  Hz 
high  frequency  cutoff.   It  was  noticed  that  a  small  rounding 
of  the  higher  amplitude  peaks  occured  in  the  waveform  with  a 
5000Hz  high  frequency  cutoff.   This  rounding  was  more  pro- 
nounced with  a  high  frequency  cutoff  of  2500  Hz.   No  apprecia- 
ble change  in  the  waveform  shape  was  observable  with  high 
frequency  cutoffs  of  20,000  Hz  and  10,000  Hz. 

Although  these  observations  are  not  conclusive  evidence 
for  neglecting  aliasing,  they  do  lead  to  the  acceptance  of  the 
calculated  power  spectra  as  being  representative  of  the  actual 
power  spectra  up  to  2500  Hz.   Furthermore,  this  analysis  indi- 
cates that  an  effort  to  increase  the  sampling  frequency  by  a 
factor  of  two  (to  10,000  samples/sec)  should  be  made.   Quali- 
tative evidence  shows  that  this  increased  sampling  rate  would 
almost  completely  eliminate  aliasing  affects  from  the  calcu- 
lated power  spectra. 

6.3   The  Fast  Fourier  Transform 

Computer  computation  of  the  amplitude  power  spectra  was 
done  through  an  intermediate  calculation  of  the  Fourier  trans- 
form coefficients.   The  relationship  between  the  .power  ampli- 
tude spectrum  and  its  associated  Fourier  transform  is  simply 
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2      2 

P   =  a    +  b 

n     n      n 

th 
where  P   is  the  n    amplitude  of  the  power  spectrum  and  a   and 
n  '  n 

b   are  the  n    Fourier  coefficients. 

n 

Instead  of  using  the  defining  equations  to  calculate  a 
and  b   as  given  in  Chapter  3,  an  algorithum  was  used  called  the 
fast  Fourier  transform  (FFT)  (Cochran,  1968).   To  use  this 
algorithum  as  supplied  in  a  Special  Scientific  Subroutine 
package  supplied  by  IBM  (IBMRef.  8)  it  is  necessary  that  the 
waveform  to  be  Fourier  transformed  consist  of  2N  samples, 
where  2N  =  2 '  and  m  is  a  positive  integer  between  3  and  20. 

The  defining  equation  for  the  FFT  is, 

2N-1 

A   =  a   +  jb   =  Z    x„exp  (-2iTjnk/N)  ;   n=0  ,  1,  ...,  N-l, 
n     n      n    ,   -  is. 

k=0 

where  A   is  the  complex  n    coefficient  of  the  FFT  and  x. 
n  k 

th 
denotes  the  k    sample  of  the  sampled  time  series. 

If  the  original  continuous  waveform  is  considered  band- 
limited  and  sampled  at  the  Nyquist  rate  it  is  possible  to 
truncate  the  waveform  after  2N  samples  and  assume  a  periodic 
continuation  of  the  2N  samples. 

The  Fourier  coefficients  of  this  series  as  computed  by 
the  FFT  differ  from  the  Fourier  coefficients  of  the  original 
bandlimited  waveform  only  by  an  aliasing  distortion. 

The  aliasing  distortion  was  shown  to  be  small  in  Sec. 
6.2  so  that  the  amplitude  power  spectra  computed  via  the  FFT 
should  be  a  realistic  approximation  of  the  actual'  unaliased 
power  spectra. 
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The  greatest  advantage  in  using  the  FFT  is  that  it  is  highly 
eff  icier.  C  in  computing  the  Fourier  coefficients  from  a  time- 
cost  standpoint.   If  the  tirr.c  series  consists  of  N  =  2   samples 
then  about  2N  log„N  arithmetic  operations  will  be  required  to 
evaluate  all  associated  FFT  coefficients.   The  number  of  oper- 
ations required  for  conventional  calculation  of  the  Fourier 

2 
coefficients  requires  at  least  N   operations.   For  example, 

in  this  thesis  the  FFT  for  N  =  8192  samples  was  computed  (i.e., 

n  =  13),  the  sum  of  the  squares  of  the  coefficients  computed, 

the  results  condensed,  normalized,  and  plotted  in  2.56  minutes. 

Typical  results  of  the  FFT  applied  to  the  nerve  and  muscle 

waveform  for  various  C0„  levels  are  shown  in  Fig.  7.1  and  Fig. 

7  2 

6.4   Cross  correlation  Function 

Two  factors  were  responsible  for  the  negation  of  this  cal- 
culation in  this  thesis.   The  first  factor  that  served  to 
eliminate  calculating  the  cross  correlation  function  between  the 
nerve  and  muscle  signals  for  various  C0_  levels  was  that  such 
a  computation  required  the  nerve  and  muscle  signals  to  be  re- 
corded simultaneously  so  that  the  CO,,  state  was  the  same  for 
both  waveforms.   The  one  channel  recorder  used  in  this  experi- 
ment did  not  allow  this. 

Furthermore,  the  calculation  of  the  cross  correlation 

2 
function  by  conventional  means  requires  at  least  N   calcula- 
tions.  With  data  for  26  experiments,  the  time  this  would  have 
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taken  would  exceed  r-  asonable  bounds.   However,  It  is  suggested 
that  a  limited  number  of  cross  correlation  functions  be  computed 
in  the  future  after  modeling  of  the  respiratory  system  has 
been  done  to  check  the  validity  of  the  model  and  a  suitable 
two  channel  FM  recorder  is  obtained. 
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CHAPTER  VI : 


ANALYSIS  OF  T:iE  NERVE  AND  MUSCLE 
POWER  AMPLITUDE  FREQUENCY  /LOTS 


7.1   Introduction 

A  respective  normalized  power  spectrum  is  depicted  in  Fig. 
7.1.   It  is  typical  of  both  the  nerve  and  muscle  power  ampli- 
tude plots  calculated  by  means  of  the  FFT  and  plotted  on  the 
IBM  360/50  computer. 

Note  that  P(f)  is  equal  to  zero  at  500  Hz  and  1000  Hz.   It 
was  necessary  to  set  the  power  at  these  frequencies  equal  to 
zero  since  unusually  large  amplitudes  appeared  at  these  two 
frequencies.   They  were  due  to  the  carrier  frequency  oscillator 
in  the  FK  tape  recorder.   The  carrier  frequency,  at  a  tape  speed 
of  .375  ips ,  was  500  Hz.   Theoretically  the  carrier  frequency 
and  any  higher  harmonics  are  eliminated  by  a  0  to  100  Hz  band- 
pass filter  built  into  the  recorder.   The  Fourier  coefficients 
of  a  sinusoid  of  frequency  f   have  a  very  large  amplitude  con- 
centrated at  f   in  the  power  amplitude  plot.   Therefore,  the 
coefficients  corresponding  to  the  FM  recorder's  carrier  frequency 
and  second  harmonic  predominated  over  the  coefficients  of  the 
waveform  at  these  frequencies. 

These  large  values  ruined  the  normalizing  routine  used  in 
the  plotting  routine.   Therefore,  they  were  arbitrarily  set  to 
zero  to  eliminate  this  problem.   Therefore,  no  physiological 
significance  is  associated  with  the  lack  of  power  at  500  and 
1000  Hz. 
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From  plots  similar  to  Fig.  7.1  dona  for  each  nerve  and 
muscle  waveform  a  cumulative  power  plot  was  made  (Fig.  7.2). 
The  purpose  of  these  being  that  the  bandwidth  of  each  signal 
could  then  be  defined. 

A  white  noise  component  was  common  to  all  the  power 
spectra  plots  as  indicated  by  the  relatively  constant  ampli- 
tudes in  Fig.  7.1  extending  from  1000  Hz  to  2500  Hz.   This  noise 
component  was  eliminated  from  the  cumulative  power  spectra  by 
a  graphical  technique.   The  bandwidth  of  a  waveform  was  defined 
as  that  frequency  below  which  90%  of  the  waveform's  power  was 
located  after  removal  of  white  noise  components.   No  frequency 
components  appeared  in  any  of  the  plots  below  30  Hz  because  the 
high  gain  amplifier  in  the  recording  system  had  a  low  frequency 
cutoff  of  35  Hz  to  eliminate  artifact  (e.g.,  body  movement  of 
probes,  heartbeat,  etc.)  from  the  recorded  xjaves  . 

Once  the  plots  were  compensated  for  the  FM  carrier  oscillator 
energy  and  white  noise  perturbations,  an  investigation  of  the  C0? 
dynamic  operator  state  and  its  affect  on  the  waveforms  bandwidth 
and  power  content  was  initiated.   The  results  of  which  are  dis- 
cussed in  the  next  two  sections. 

7.2   The  Effect  of  C02  Dynamic  Operator  State  on  the  Muscle 
Waveform  Power  Spectra  Characteristics 

The  two  most  apparent  affects  for  various  operator  states 
were  found  to  be  the  waveform's  bandwidth  and  its  relative  power 
content.   The  relative  power  content  was  based  on  the  amplitude 
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of  the  cumulative  power  spectra  at  Its  horizontal  asymtote 
(after  accounting  for  normalization).   The  largest  asymptote 
was  defined  as  100%. 

The  state  of  the  dynamic  operator  was  indicated  by  the 
measure  of  inspired  C0?.   The  results  of  bandwidth  analysis  are 
shown  in  Fig.  7.3.   Fig.  7.3  graphically  illustrates  how  the 
C0„  level  affected  the  muscle  waveform  bandwidths  in  the  experi- 
ment.  No  bandwidth  values  appear  below  6%  C0„  because  all 
muscle  action  ceased  and  the  associated  electrical  activity  was 
not  great  enough  to  be  distinguished  from  noise  in  the  recording 
sys  tem. 

Fig.  7.3  indicated  that  decreasing  levels  of  C02  still  main- 
taining phasic  respiratory  movement,  produce  an  almost  linear  de- 
crease in  muscle  waveform  bandwidth.   In  Fig.  7.3,  changing  the 
state  of  the  operator  from  14%  CO   to  10%  CO   caused  a  corres- 
ponding decrease  of  bandwidth  from  194  Hz  to  153  Hz. 

At  a  level  of  8%  C0„,  the  animal  was  observed  to  be 
approaching  apnea.   When  this  occurred,  Fig.  7.3  indicated  a 
comparitively  large  bandwidth  (217  Hz).   At  6%  C0?  the  muscle 
bandwidth  decreased  to  177  Hz  and  below  6%  C0_  no  muscle  wave- 
form was  observable. 

Fig.  7.4  is  a  plot  of  the  relative  power  of  each  muscle 
waveform  versus  C0„  level.  It  is  observed  that  the  relative 
power  decreased  as  the  C0„  level  decreased. 
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7.3   The  Effect  of  CO-  Dynamic  Operator  State  on  the  Nerve 
Waveform  Power  Spectra  Characteristics 

Results  of  the  adjusted  cumulative  power  spectra  plots  are 
summarized  in  Fig.  7.5  for  the  nerve  waveform.   The  bandwidth  is 
plotted  as  a  function  of  the  inspired  CO-  level.   Decreasing  C0_ 
levels  have  an  entirely  different  affect  on  the  nerve  waveform 
bandwidth  than  they  did  on  the  muscle  waveform  bandwidth  (at 
least  during  phasic  respiration). 

Notably,  as  the  C0„  level  decreased  from  14%  to  10%  the 
nerve  bandwidth  remained  practically  constant.   The  bandwidths 
were  311  Hz,  321  Hz,  and  316  Hz  corresponding  to  CO-  levels  of 
14%,  12%  and  10%  respectively. 

Fig.  7.4  indicates  that  as  the  C02  level  decreased  from  14% 
to  10%,  the  relative  power  in  the  nerve  waveform  displayed  a 
rising  trend.   The  rising  trend  is  diametrically  opposite  to  the 
decrease  in  relative  power  displayed  by  the  muscle  waveform  for 
the  same  CO-  levels. 

Again,  as  a  state  of  apnea  was  approached  (8%  CO-)  a 
dramatic  change  in  the  power  characteristics  of  the  nerve  wave- 
form transpired.   Abruptly,  both  the  bandwidth  and  relative  power 
increased  to  near  maximum  values  at  the  onset  of  apnea.   The 
relative  power  reached  a  maximum  at  6%  C09.   Maximum  bandwidth 
(817  Hz)  occurred  during  apnea  at  a  CO-  level  of  4%. 

C09  levels  of  2%  and  0%  displayed  a  decreasing  trend  in 
both  the  nerve  waveform  bandwidth  and  relative  power.   Notice 
that  at  low  levels  of  CO-  during  which  no  respiratory  muscle 
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action,  existed  there  was  still  quice  a  bit  of  electrical 
activity  in  the  nerve.   In  fact,  Fig.  7.4  shows  that  the 
relative  power  at  the  lower  C0„  levels  is  as  great  or  greater 
than  the  relative  power  associated  with  higher  levels  of  C02 
producing  phasic  respiration. 

7 .  4   Conclusion 


The  discussion  of  the  muscle  waveform  power  characteristics 
in  Sec.  7.2  suggests  that  at  least  two  mechanisms,  or  two 
different  dynamic  operators  are  associated  with  the  respiratory 
system  response  to  various  C0„  states.   The  change  from  one 
operator  to  the  other  occurring  at  that  level  of  C0„  at  which 
phasic  respiration  ceases  and  apnea  begins.   It  also  appears 
likely  that  the  affect  of  the  dynamic  operators  state  on  the 
muscle  waveform  during  phasic  respiration  is  to  produce  a  lin- 
early related  muscle  waveform  bandwidth  and  relative  power. 

The  abrupt  changes  in  the  nerve  waveform  power  characteristics 
discussed  in  Sec.  7.3  also  imply  the  existence  of  two  mechanisms, 
or  two  different  dynamic  operators  associated  with  C02  dynamic 
operators  state  and  its  affect  on  the  respiratory  nerve  elec- 
trical activity.   For  phasic  respiration,  the  most  significant 
factor  the  state  of  the  dynamic  operator  has  on  the  nerve  wave- 
form is  its  relative  power  content.   The  relative  power  content 
increasing  as  the  C09  level  decreases. 


53 


7.5   Recommendations  for  Further  Research  in  this  Area 


The  next  step  in  the  continuing  study  of  nerve  and  muscle 
behavior  in  the  respiratory  system  should  be  done  with  the  cor- 
relation functions  defined  in  Chapter  V.   The  results  of  which, 
through  the  analysis  of  Sec.  5.8,  can  lead  to  a  mathematical 
model  of  the  respiratory  system  using  the  model  derived  in 
Chapter  III.   The  first  model  should  try  to  account  for  the  ob- 
served respiratory  amplitude  changes  produced  by  various  states 
of  the  dynamic  operators.   After  this  is  done  a  model  producing 
waveforms  similar  in  bandwidth  and  relative  power  to  those  in  this 
thesis  might  be  produced  by  a  digital-analog  computer  facility. 

The  apparent  relationship  between  the  state  of  the  C02  dy- 
namic operator  and  the  nerve-muscle  waveform  descriptors  mentioned 
previously  leads  to  an  exciting  idea  suitable  for  further  in- 
vestigation.  Physiologists  have  been  seeking  a  way  to  characterize 
electrical  activity  of  the  muscle  under  tension  according  to  the 
amplitude,  shape,  and  frequency  of  occurance  of  individual  muscle 
cell  potentials.   The  integrated  result  of  many  muscle  cells 
firing  under  controlled  levels  of  tension  is  identical  to  the 
muscle  waveform  recorded  in  this  experiment. 

Therefore,  since  the  state  of  the  CC-   dynamic  operator 
controls  muscle  tension,  and  waveforms  similar  to  those  analysed 
in  this  experiment  represent  the  combined  effect  of  the  indi- 
vidual potentials  (whose  character  is  to  be  described),  can  the 
spectra  descriptors  i.e.,  bandwidth  and  relative  power  be  used 
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for  this  purpose?   This  is  an  intriguing  question  and  when 
answered  would  greatly  increase  the  knowledge  of  muscle  action. 

Hardware  for  the  digital-analog  system  interface  was  part 
of  the  effort  of  this  thesis.   Therefore,  research  in  the  future 
can  be  directed  to  higher  theoretical  levels  and  more  sophisticated 
models  of  the  respiratory  system.   Above  all,  it  is  highly  recom- 
mended that  many  more  experiments  be  performed  to  confirm  the 
conclusions  of  this  thesis  and  to  open  avenues  to  the  study  of 
the  respiratory  system  that  are  not  dreamed  of  now. 
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APPENDIX  A 
ANALOG-TO-DIGITAL  CONVERSION  SYSTEM 

A.l   Introduction 

This  appendix  includes  the  majority  of  work  done  on  the 
design  and  development  of  the  analog-to-digital  conversion 
system.   The  data  processing  system  represents  the  majority 
of  work,  time,  and  effort  in  the  realization  of  this  thesis. 
At  the  time  of  its  conception  a  need  for  analog-to-digital 
conversion  was  desperately  needed  by  several  departments 
at  Kansas  State  University.   Addition  of  an  IBM  360/50 
at  the  university  computing  center  opened  the  door  to 
extensive  and  sophisticated  analysis  of  experimental  data. 
Unfortunately,  almost  all  experimental  data  was  analog 
in  nature  and  not  compatible  with  digital  computer  input 
processes  (e.g.,  nerve  and  muscle  signals). 

Another  area  where  a  data  conversion  system  was  needed 
was  in  the  Department  of  Psychology.   They  were  involved 
in  a  human  tracking  experiment.   It  was  found  that  hand 
scoring  of  raw  analog  data  varied  widely  according  to  the 
individual  scoring  a  particular  experiment.   Besides  this, 
handscoring  was  a  laborious  task  involving  at  least  one 
afternoon  x.o    score  only  one  trial  of  an  experiment. 
Obviously,  something  had  to  be  done.   The  data  processing 
system  described  herein  was,  in  part,  a  solution  to  their 
problems.   For  further  information  on  this  subject  see 
Dale  E.  Bentrup's  Master's  thesis,  Kansas  State  University, 
June  1,  1968.   Using  this  system  Bentrup  has  reduced  the 
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scoring  time  to  a  minimum  and  greatly  enhanced  the 
reliability  of  the  experimental  results.   In  fact,  he 
can  now  score  the  same  number  of  trials  in  one  afternoon 
that  would  have  taken  a  month  to  score  by  hand. 

This  writer's  interest  was  in  the  application  of 
information  theory  and  statistics  to  biological  data 
analysis.   Data  was  available  to  work  with,  but  it  was 
in  analog  form,  and  not  suitable  for  computor  methods 
of  analysis.   The  following  sections  therefore  include 
the  design  and  development  of  a  data  processing  system 
that  enabled  the  theoretical  work  to  commence. 

A. 2   Analog-to-Digital  Converter  Output  Number  System 

The  combination  multiplexer  and  analog-to-digital 
converter  (ADC)  used  in  the  system  is  a  Series  QkB    ADC 
built  and  purchased  from  Texas  Instruments.   It  operates 
in  straight  binary  format.   It's  output  format  is  either 
offset  binary  or  twos  complement. 

An  offset  binary  system  is  defined  as, 
+  full  scale  =  1  1111  1111  1111  1 
zero  =  1  0000  0000  0000  0 

-  full  scale  =  0  0000  0000  0000  0 

and  a  twos  complement  number  system  is  defined  as  , 
+  full  scale  =  0  1111  1111  1111  1 
zero  =  0  0000  0000  0000  0 

-  full  scale  =  1  0000  0000  0000  0 

Notice  the  only  difference  between  the  offset  binary 
and  twos  complement  number  system  is  that  the  sense  of  the 


sign  bit  is  reversed  in  the  twos  complement  system.  The 
offset  binary  system  option  was  chosen  to  facilitate  the 
reconstruction  of  the  sampled  voltage  in  the  IBM  306/50. 

For  the  offset  binary  system  the  input  voltage  can 
be  reconstructed  from  the  binary  output  bits  with  the 
formula , 


e  . 

i 


(A.l) 


where  ,    e.       =    input    voltage 

E  =    machine    full-scale    voltage    10V 

n         =    number    of    output    bits    13 

A     -•  /  \ 

k    =  binary  output  bits  (A   =  sign  bit;. 

The  resolution  of  the  data  conversion  system  is  equal 
to  the  minimum  Voltage  level  discernable  in  the  binary 
output.   In  this  case,  the  minimum  voltage  discernable 
corresponds  to  the  value  of  the  LSB,  where, 

E 

LSB  =  -f|  (A. 2) 

E,   is  10  volts  in  the  Series  8U8  ADC  so  that, 
fs 

LSB  =  .00121+7  V  (A. 3) 

The  accuracy  of  voltage  resolution  is  approximately 

one  thousandth  of  a  volt  as  can  be  seen  in  Eq.  A. 3- 
A. 2.1   Multiplexer  Operation 
Incorporated  in  the  same  package  as  the  ABC  is  a 

multiplexer  capable  of  switching  up  to  l6  analog  signal 
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inputs  to  a  single  output.   The  switching  sequence  is 
selectable  and  may  he  i.i  either  of  two  normal  operating 
modes,  sequential  or  addressable.   In  tie  sequential  mode, 
input  channels  are  sampled  by  internal  clock  circuitry  or 
are  sampled  at  a  rate  determined  by  external  equipment. 
In  the  addressable  mode,  input  channels  may  be  sampled 
in  any  desired  sequence  controlled  by  a  seven-bit  channel 
address  from  external  equipment  (Texas  Instruments, 
pp.  1-3). 

In  the  addressable  mode,  an  address  can  be  entered 
only  when  the  front  panel  switch  "mode"  is  in  the  address- 
able position  and  the  "clock"  switch  is  in  EXT  (external). 
This  places  the  instrument  in  addressable  mode  and  turns 
timing  control  over  to  external  equipment.   An  address 
is  entered  by  providing  the  desired  binary  address  at  the 
proper  input  pins  and  then  providing  a  -6V  pulse  of  at 
least  two-microseconds  duration  at  the  address  set  input 
pins  (Fig.  A.k).       When  the  address  has  been  set  (1.5  psec. 
after  receiving  an  address  set  command)  an  internally 
generated  -6V  pulse  is  generated  and  gated  to  the  clock 
output  pin  (Fig.  A.U). 

A. 2. 2   Analog-to-Digital  Conversion 

The  clock  pulse  occurring  after  the  address  has  been 
set  is  connected  to  the  "sample  initiate"  command  input 
pin  (Fig.  A.M.   This  clock  pulse  begins  the  conversion 
process.   The  clock  pulse  closes  an  input  gate  applying 
the  input  signal  to  a  holding  capacitor  through  a  buffer 
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stage.   The  input  signal  is  sampled  for  11  ysec.  and  then 
the  input  gate  opens  leaving  the  input  signal  trapped  on 
the  holding  capacitor. 

The  held  signal  is  applied  as  one  input  to  a  differ- 
ential comparator.   The  second  input  is  taken  from  an 
internal  precision  voltage  supply.   The  internal  supply 
produces  an  output  voltage  directly  proportional  to  the 
"binary  number  contained  on  the  output  pins.   Thus,  if  the 
contents  of  the  binary  output  pins  are  manipulated  so 
the  two  inputs  to  the  differential  comparator  are  equal 
(i.e.,  the  comparator  output  is  zero),  then  the  binary 
number  on  the  output  pins  is  directly  proportional  to  the 
input  signal  stored  on  the  holding  capacitor.   The  tech- 
nique used  to  manipulate  the  output  binary  bits  to  represent 
the  held  voltage  is  called  the  successive  approximation 
technique.   An  excellent  discussion  of  this  technique  is 
found  on  page  3-13  of  The  Texas  Instruments  instruction 
Manual  No.  176766-OOOI,  Series  81; 8  ADC  Manual.   Thirty- 
nine  usee,  are  required  for  the  successive  approximation 
of  the  held  signal.   At  the  end  of  this  period  the  output 
bits  are  available  to  external  equipment  (Fig.  A. 3)  as  a 
parallel  ik    bit  display. 

Basically,  four  events  must  occur  for  conversion: 

1.  Analog  signals  connected  to  multiplexer  inputs 

2.  Channel  address  entered  in  multiplexer 

3.  Address  set  command  pulse  transmitted  to  multiplexer 
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k.       Clock  pulse  routed  to  sample  initiate  command  input. 

The  binary  output  appears  39  usee,  after  initial  address 

set  command. 

A. 2. 3   Electrical  Specifications  of  ADC  Control 

Signals 

The  specifications  are  tabulated  in  table  A.l  below. 

Multiplexer 

Number  of  input  channels       l6 

Signal  Input  -  10  volts 

Logic  levels  1  =  6V,  0=0  Volts 

Channel  Address  seven-line,  single  rail, 

binary  address 

Address  set  command  external  clock 

Clock  (internal)  -6V  pulse  occurring 

1.5  usee,  after  address 
set  command 

Input  impedence  -  lOO.Mftshunted  by  -  25pf 

Maximum  source  impedence  1000  ft 

Analog-to-Digital  Convertor 

Total  system  cycle  time  39  usee. 

Parallel  digital  output  single  rail,  Ik    bits 

Control  Signal  Requirements 

< 

Rise  time  -  .5  usee. 

Pulse  duration  -  2  usee. 

< 

Fall  time  -  o  usee. 

Maximum  current  +   2.0  ma 

TABLE  A.l.  Electrical  specifications  of  ADC 
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A. 3   Digital  Stepping  Recorder  (D3R) 

The  digital  stepping  recorder  used  in  the  data  processing 
system  is  a  Precision  Instruments  DSR  1^20.   Its  general 
specifications  are  as  listed  in  table  A. 2.   Electrical 
specifications  are  listed  in  table  A. 3.  (Digi.  Data  Corp.) 

1.  Recording  Rate  200  characters /sec . 

2.  Interrecord  Gap  <130  ms 

3.  Tape  Width  3-S  inch  standard 

h.       Number  of  Tracks  7 

5.   Interrecord  Gap  Generation  3M  inch  (under  ext. 

control) 


6 

8 

9 


Bit  Density 
Parity 


200  BPI 

longitudinal  and  lateral 
selectable  odd  or  even 


End  of  File  Gap  and  Mark   3.^  inch  (controlled 

from  front  panel) 

Tape  Load  Point  Positioning   Photocell 

TABLE  A. 2   General  specifications  of  DSR 


3 
h, 


Signal  Inputs 
Data  Inputs 

Step-Record  Command 


Interrecord  Gap  Command 


Binary  1  =  >-6V, 
Binary  0  =  0V 

7  lines  ,  must  be  present 
during  stop-record  cmd. 
pulse;  7-5K  impedance 

50-100  psec.  negative 
pulse;  7-5K  impedance 

20  ps  negative  pulse  , 
15K  impedance 


TABLE  A. 3   Electrical  Specifications  DSR 
6.3.1   Operation  of  the  DSR  1^20 
Operation  is  very  simple.   3ev:.-n  lines  of  parallel  data 
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are  brought  to  the  recorder  input  (Fig. A. 9).   A  negative 
pulse  of  approximately  80  psec.  is  applied  to  the  3tep  record 
input.   This  automatically  writes  the  input  and  advances 
the  tape.   The  one  limitation  that  can  not  he  exceeded  is 
the  synchronous  digital  stepping  rate  i.e.,  200  steps/sec. 
Exceeding  this  by  a  large  factor  would  be  disasterous  to 
the  stepping  mechanism. 

k.h      Wiring  Diagrams  Rack  Mounted  Analog-to-Digital 

Converter  and  Digital-to-Analog  Converter 
Tc  facilitate  operation,  the  ADC  unit  was  mounted  in 
an  upright  rack  (Fig.A.l)  along  with  a  companion  Digital- 
to-Analog-Converter  DAC.   The  DAC  was  not  used  in  the  data 
processing  system  described  in  this  thesis.   It  is  antici- 
pated that  it  will  be  useful  at  a  later  date  as  a  link 
between  the  ATHENA  digital  computer  and  an  analog  computer. 
Therefore,  this  section  will  include  all  wiring  documentation 
of  the  rack  mounted  ADC  and  DAC. 

All  electrical  signals  on  unit  chassis  connectors 
J117,  J118,  and  J119C  of  the  ADC  and  those  on  Jllo  of  the 
DAC  were  brought  to  Junction  Box  2.   J-Box  2  is  located  at 
the  boitoa  of  the  rack  and  easily  accesible  to  underfloor 
wiring  and  cabling  (Fig.  A. 2).   Both  units  can  be  com- 
pletely controlled  by  external  equipment. 

A  number  of  test  and  monitor  signals  were  brought 
to  special  panels  on  the  front  of  the  rack  from  the  ADC 
and  DAC.   Refer  to  Figs.  A.3-A.10  for  specific  definition 
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of  these  signals.   Note  that  only  two  ADC  analog  input 
channels  (zero  and  one)  ere  brought  to  the  upper  front 
panel  in  Fig.  A.l.   These  are  especially  useful  for  dig- 
itazation  from  external  sources  with  only  one  or  two 
channels.   Otherwise,  the  analog  inputs  on  J36  (Fig. A. 2) 
must  be  used. 

A. 5   The  Digital  Computer  and  Program  for  Analog- 
To-Digital  Tape  Conversion 

The  computer  used  to  control  all  peripheral  data 
conversion  equipment  was  a  USAF  Model  SM  68  Missile  Weapon 
System  semi-general  purpose  digital  computer.   It  was  built 
by  Remington  Rand  Univac  for  use  in  the  Titan  missile  base 
at  Lowry  Air  Force  Base  in  Denver,  Colorado. 

When  the  site  was  deactivated  by  the  Air  Force  the 
computer  was  obtained  by  the  Department  of  Electrical 
Engineering  at  Kansas  State  University.   At  the  University 
the  SM  68  Missile  computer  is  called  the  Athena.   Athena 
has  a  relatively  small  memory  of  256  twenty-four  bit  words. 
A  magnetic  drum  storage  section  is  available  which  may  store 
up  to  8192  seventeen  bit  words.   The  drum  storage  provides 
permanent  storage  for  program  instruction  and  constants. 

Once  a  program  is  started  it  proceeds  sequentially 
through  the  8192  addresses,  excuting  each  instruction 
before  preceding  to  the  next  address.   The  minimum  time  for 
one  instruction  to  be  excuted  is  kO    psec.   For  more 
information  about  the  programming  and  operation  of  Athena 
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refer  to  the  AF  Technical  Manual  A?Ck(6h5)-20    and  AFOfc 
(69*0-71  titled,  SM  63  Missile  Weapon  System  Rail?.- *5£.£.1L?L 
5.H.L^S.£.^.  S^s_t_e_m  Co  input  or .   The  following  sections  will 
assume  a  knowledge  of  sections  V,  VI,  VII,  VIII,  IX,  and 
X  of  the  technical  manual.   These  sections  include  computer 
control,  arithmetic,  magnetic  core  storage,  magnetic 
drum  storage,  programming  instructions,  and  input-output 
data  routing.   All  instructions  in  the  following  sections 
are  given  in  octal  code  as  are  all  drum  addresses  and 
constants  stored  in  core  locations. 

A. 5.1   Input-Output  Signals  and  Bit  Positions  On  Control 
Console 

The  data  conversion  system  requires  four  command 
signals  from  the  Athena  computer.   They  consist  of: 

1.  a  multiplexer  address  in  straight  binary  code 

2.  an  address  set  command  for  the  multiplixer 

3.  a  step-record  command  for  the  digital  tape  deck 

h .   an  end  of  record  command  for  the  digital  tape  deck. 
In  addition,  Athena  must  act  as  a  buffer  for  the  1*;  binary 
bit  ADC  output.   Figs.  A. 3  -  A. 10  are  wiring  diagrams 
showing  the  connections  of  the  peripheral  equipment  with 
the  Athena  computer.   Fig.  A. 11  is  a  listing  of  the  computer 
registers  as  they  appear  on  the  control  console  with  the 
corresponding  four  command  signals  and  Ik    bit  binary 
sample  inputs.   Fi^s.  A. 12  and  A. 13  display  the  wiring 
cables  at  the  Athena  max-.ng  connector  and  J-3ox  8030C. 
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117-  1"  • 

Sign 

2-  •  • 

Sign 

3-  •  • 

Bit  1 

1|  •  •  • 

Bit  1 

5... 

Bit  2 

6  •  •  • 

Bit  2 

?••  • 

Bit  3 

8  •  •  • 

Bit  3 

Q  .  .  . 

Bit  4 

10  •  •  • 

Bit  4 

11"  • 

Bit  5 

12 


3it    5 


13* • •    Bit    6 


•J30-1 

■J36-2 

■NC 

•J36-5 

•NC 

■J36-7 

•NC 

■J  36 -9 

•NC 

•J36-11 

•NC 

•J36-13 
•NC 


14»««    Bit    o     •••J3o-15 


15 


Bit    7 


lo  •  •  • 

Bit 

7 

17"  • 

Bit 

B 

18  — 

Bit 

8 

19*  •  • 

Bit 

9 

20'  •  • 

Bit 

9 

2i.  .  . 

Bit 

10 

22 
23 


Bit    10 
Bit    11 


•NC 

•J36-17 

•NC 

•J36-19 

♦NC 

•J36-21 

•NC 

*J36-23 

•NC 


P117-24 Bit  11 

•  •  •  J36-25 

25  •  •  •  Bit  12 

•  •  -NC 

26  •  •  •  Bit  12 

• ..J36-27 

27 Bit  13 

•  •  «NC 

28-  •  •  Bit  13 

• • -J36-29 

29  •  •  •  Bit  14 

•  •  'NC 

30 • • •  Bit  14 

•  •  «J36-31 

31« • •    NC 

•  •  •  N  C 

32 • • •    NC 

•  •  «NC 

33- ••    NC 

•  •  «NC 

34- ••    NC 

•  •  «NC 

,c     Digit 

Complete 

• • «TB1-1 

36--.  Digit 
J     Complete 

•  •  •  H  C 

Pulse 
i{           Width 

T31-2 

38...  Pulse 
*             Width 

•  •  -NC 

Anerture 
^         Control 

•  •  -T31-3 

40  •  •  •    NC 

•  •  -:;c 

41«  ••    NC 

•  ••NC 

l+2«  •♦    NC 

•  •  •  N  C 

43 • • «AD  Low 

•  •  -TB1-4 

44  •  • 'AD  High 

•  •  •T31-5 

45 • • •    NC 

•  •  -NC 

46 • • • Ground 

•  •  -Sheilds 

Fig. A. 3.       Wiring    assignments    for    connector   P117    on    ADC 


P118-    l-« 

2-  < 

3" 
4., 

5" 
6- 

7- 

8- 
9- 

lO- 
ll- 
12- 
13- 
14- 


NC 

„6 


2° 

NC 

„6 


•  •  •  N  C 

• •  -J35-1 
•• -J35-3 

•  •  «J35-5 
• • -J35-7 

•  •  -J35-9 

•  •  -J35-H 

•  -  -J35-13 

•  •  •  N  C 

• • -J35-15 
• • -J35-17 
• • * J35-19 

•  •  -J35-21 

•  --J55-23 


P118-24 


25 


26.  • 


27- 

28- 

29- 

30- 

31« 

32- 

33' 
34' 

35' 

36 

37 


Frame 
'Output 

•  •  .TB4-11 

Mode  Se2 
Cmnd 

■ect.T34-4 

•   NC 

•  •  «NC 

•   NC 

•  •  •  N  C 

•   NC 

•  •  'NC 

•   NC 

•  •  -NC 

•   NC 

.  •  • 'NC 

•   NC 

•  •  'NC 

•   NC 

•  •  •  N  C 

•   NC 

•  •  'NC 

•   NC 

•  •  -NC 

•   NC 

•  •  'NC 

•   NC 

•  •  'NC 

•   NC 

•  •  -NC 

(Continued  Fig. A.1*.) 
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15 

16 

17 
18 
19 
20 
21 
22 
23 

Fig. A 


Addres  s 
Set  Rdy 
.•ess 
Set  Cmnd 

Sample 
In.  Rdy 
Sample 
In .  Cmnd 
Mode  Cnfrm. 

Output 
Ext .  Reset 
Cmnd 
Clock 

*  • 

Output 


•J35-25 
•J35-27 
•TBU-7 


TBl-l 
TBU-8 
TB**-2 
:TBI-9 
TB^-3 
T3H-10 


33-- 

NC 

•  -NC 

39" 

• 

•  -NC 

UO  •  • 

NC 

•  -N  C 

in.. 

u  c   ' 

•  'NC 

k2- 

NC 

•  «NC 

k  3  •  • 

NC 

•  «NC 

Ul*  -  - 

NC 

•  -NC 

U5«- 

!IC 

•  •  »NC 

46  •  • 

Logi  c 
Ground 

•  •  -T34-6 

k.      Wiring  assignments  for  connector  Pll8  on  ADC 


P119C-  1 Channel 

2 

3 

k 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 Common 

18 Common 

19 Common 


0  TB3-10 

1  TB3-11 

2  J37-5 

3  J37-7 

h    J37-9 

5  J37-11 

6  J37-13 

7  J37-15 

8  J37-17 

9  J37-19 

10  J37-21 

11  J37-23 

12  J37-25 

13  J37-27 

Ik    J37-29 

15  J37-29 

TB3-9 

TB3-9 

J37-33 


Fig. A. 5.   Wiring  assignments  for  connector  P119C  on  ADC 
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P117-35- 

TB1-1  • 

•J36-33&TB2-1 

TB1-1  • 

TB2-1  • 

•J31 

37- 

2  • 

•J36-35&TB2-2 

TB1-2  • 

2  • 

•J32 

39- 

3  • 

•J36-37&TB3-1 

T31+-7  ' 

3  • 

•J26 

1+3- 

I    - 

•J36-39&TB3-2 

TBfc-8  • 

1+  - 

•J2  7 

kk- 

5  • 

•J36-ltl&TB3-3 

TBlt-9  • 

5  • 

•J26 

hb- 

o     • 

•J36-U3 

NC 

6  - 

•  NC 

NC- 

7  • 

NC 

TBlt-10  • 

7  • 

•J29 

NC- 

3  - 

•   NC 

TB>*-11- 

8  • 

•J30 

NC- 

9  • 

•   NC 

NC 

9  • 

•  NC 

NC- 

10- 

•   NC 

NC 

lO- 

•  NC 

NC-  - 

11- 

-   NC 

NC 

ll- 

•  NC 

TB1-3   • 

TB3-1  • 

•J23 

P118-18 

•  -TBU-l- 

•J35-29&TB3-^ 

TB1-U   • 

2  • 

•J2fc 

P118-20 

2- 

•J35-31&TB3-5 

TB1-5 

3  • 

•J25 

P118-22 

3- 

•J35-33&TB3-7 

TBU-1 

1+  • 

•J19 

P118-25 

k. 

•J35-33&T33-8 

TBH-2 

5  • 

•  J20 

NC 

5. 

•  NC 

NC 

6  - 

•  NC 

P118-U6 

r 

O  ♦ 

•J35-37 

TB1+-3 

7  • 

P118-17 

7. 

•J35-39&TB2-3 

TBk-k 

8  - 

•J22 

P118-19 

8- 

•J35-H0S:T32-i+ 

P119C-17 

9  • 

•P119C-13 

P118-21 

9. 

•J35-1+lS:T32-5 

P119C-1- 

lo- 

•J37-1&J33 

P118-23 

10 

.J35-1+2&TB2-7 

P119C-2- 

ll- 

•J37-3&J3^ 

P118-2U 

11 

•J35-i+3&TB2-8 

Fig. A. 5.   Wiring  assignment  for  terminal  boards  in  J-3ox  1 


J35-1 


i 
k 


7 

3 

9 
1C 

11 
12 

13 


0 


NC 
25 
NC 
2k 

NC 

•5 
2 

NC 

22 

NC 

21 
NC 

.-0 


P118-2 

NC 
P118-3 

NC 

P118-U 

NC 

P118-5 
NC 

P118-6 
NC 

P118-7 
NC 

P118-8 


J35-21+  • 

-  •   NC   •  •  • 

NC 

25  • 

..   21   ... 

P118-15 

26  - 

•  •   NC 

NC 

27  • 

..   2o   ... 

P118-16 

28  • 

•  •   NC 

NC 

29  • 

Add. Set 

>  *              •  •  ♦ 

Cmnd 

TBU-1 

30  - 

NC 

NC 

31  •• 

Sample 
In  C mn  d  . 

TB1+-2 

32  • 

,  .   NC 

NC 

33  • 

^Ext .Reset 
Cmnd 

T3U-3 

31+  • 

,  .   NC 

NC 

35  • 

Mode  Set. 

•               •  •  • 

Cmnd 

TBfc-U 

36  - 

•   NC 

NC 

(Continued  Fig. A. 6.) 


]  5 

lb 

17 

18 

19 

20 

21 
22 
23 


NC 
2  5 

NC 
2k 

NC 

23 
NC 
22 

i  g  .  A  .  6 
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NC 

P118-10 
NC 

P118-11 

NC 
P118-12 

NC 

P118-13 

NC 
P118-14 


37 

38 
39 

1*0 

1*1 

42 
1*3 

45 


Logi  c 
Groun d 

KG 
Add. Set 

Re  ady 

Gamp  e 
In.l 

Mode 
Con  fi  rm 

Clock 
Output 

Frame 
Output 

NC 

NC 


TB4-6 

NC 
TBI*  -7 

TBI* -8 

T  B  It  -  9 
TBI* -10 

TB4-11 

IIC 
NC 


Wiring    assignments    for    J35    on    J-3ox    2 


J36-1 

2 
3 
4 

5 
6 

7 
8 
9 

10 
11 

12 

13 

14 

15 
16 
17 
18 
19 
20 
21 
22 
23 


Sign 

NC 
Sign 

NC 
Bi-c  1 

NC 
Bit  2 

NC 
Bit  3 

NC 
Bit  1* 

hi  o 


Bit  5 
NC 

Bit  6 

NC 
Bit  7 

NC 
Bit  8 

NC 
Bi-c  9 

NC 
Bit  10 


P11T-1 

NC 
P117-2 

NC 
P117-3 

NC 
P117-5 

tt  n 

P117-7 

NC 
?ll'7-9 

NC 

P117-H 

NC 

P117-13 

NC 
P117-15 
NC 

P117-17 

NC 

P117-19 

NC 
P117-21 


21*  .. 

NC   ... 

NC 

25  .. 

.  Bit  11.  . 

P117-23 

26  .  . 

.   I.  C   ... 

NC 

27  •  • 

.  Bit  12  .  . 

P117-25 

28  .  . 

N  C   ... 

NC 

29  •  • 

.  Bit  13. • 

P117-27 

30  .  . 

.   NC 

NC 

31  .  . 

.  Bit  14.. 

.P117-29 

32  .  . 

.   NC 

NC 

33  .  . 

.Dig. Compl 

.TB1-1 

31*  .. 

N  C   ... 

NC 

35  .  . 

Sample         2 
'Pulse  Width 

36  .. 

.   NC 

NC 

37  •  . 

Aperture 
'  Cntr  ••' 

TB1-3 

38  .. 

NC 

NC 

39  .  . 

.AD  Low . . . 

TB1-1* 

ho    .  . 

.   NC 

NC 

1*1  .. 

.AD  High. . 

T31-5 

H2  .  . 

NC   ... 

NC 

43  .. 

. Ground . . . 

TB1-6 

1*1*  .  . 

.   NC 

NC 

45  •• 

.   NC 

NC 

Fig. A. 7.   Wiring  assignments  for  connector  J36  on  J-.oox  2 
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J37-1 TB3-10 

2 N  C 

3 TB3-11 

4 NC 

5 P119C-3 

6 PI19C-6 

7 P119C-7 

Note:   The  regaining  pins  ( 8  —  i+  5  )  on  J37  are  not  connected 


Fig. A. 8.   :Wiring  Assignments  for  connector  J37  on  J-Eox2 


J38-1 
2 
3 
I* 

5 
6 

7 
8 
9 
10 
11 
12 
13 
14 

15 
16 

17 
18 

19 
20 

21 

22 


Bit 

0  . 

. 

P116-1 

Bit 

*} 

. 

Pllo-2 

Bit 

2  . 

. 

Pllo-3 

Bit 

3  • 

. 

PllO-4 

Bit 

4  . 

. 

P116-5 

Bit 

5  • 

. 

P116-6 

Bit 

6  . 

. 

P116-7 

Bit 

7  • 

* 

P116-8 

Bit 

3  . 

a 

P116-9 

Bit 

9  . 

, 

P116-10 

Bit 

10. 

. 

P116-11 

"D  J  — 

11. 

# 

P116-12 

Bit 

12. 

m 

P116-13 

Bit 

13. 

. 

P116-14 

NC 

. 

. 

NC 

Common . 

« 

P116-16 

Common . 

. 

P116-17 

Hi  xt 

.Ref 

■ 

Pll6-l8 

Ext 

.Ref 

Gnd.  Pll6 

Reg 

.Set 

P116-20 

Ch 

Reg, 

Ch, 

Reg, 

Ch, 


1 
Set 
2  ' 
Set 


■19 


.  P116-21 
.  P116-22 


24  ., 

.  . .  J  ... 

NC 

25  . 

.  N  C  ... 

NC 

26  . 

.  NC  ... 

NC 

27  • 

•    -  1  -j         •  •  • 

NC 

28  . 

.  NC  ... 

NC 

29  • 

.  NC  ... 

NC 

30  . 

.  NC  ... 

NC 

31  . 

Ch.lGND. 

P117&J1 

32  . 

Ch . 1  Out 

P117&J1 

33  . 

Ch.2GND. 

P118&J2 

34  . 

Ch.2  Out 

P118&J2 

35  . 

Ch.3GND. 

P1192:J3 

36  . 

Ch.3  Out 

P119iJ3 

37  • 

.  NC  ... 

NC 

38  . 

Ch .4Gnd. 

P120&J4 

39  • 

Ch . 4  Out 

P120&J4 

40   . 

Ch .5 Gnd 

P121&J5 

41  . 

Ch.5  Out 

P121&J5 

42  . 

Ch .6Gnd 

P122&J6 

43  . 

Ch.6  Out 

P122&J6 

44  . 

.  NC  ... 

NC 

Note:Jl-Jo  are  on  the  front  panel 


Fig. A. 9 


Wiring  Assignments  for  Pll6  on  the  BAC  and  J38  on  J-Box2 


(Continued  Fig. A. 9.) 
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Athena 

Plug-Pin 

Number 


P10105-1+1 
1+2 
1+3 
1+1* 
hi 
1*6 
hi 
1*8 

k9 
50 

51 
52 
53 
5h 
55 
56 

57 
P10205-1 
2 
3 
h 
5 
6 

7 

£ 

9 

10 

11 

12 

13 

Ik 

15 

16 

17 
18 
19 

2C 
21 
22 
23 
24 
25 
26 


Athena         .80  300 
Logic      Terminal  Board 
Name         Assignments 


27 
28 

29 


Ace . R-06 

Acc.R-07 
Acc.R-08 

Acc.R-09 
None 
.None 

.Targ.Ident .C 

.  Common 

,Targ. Ident .C 

.Common 

.Miss . Distanc 

,  Common 

.Plotting  3D 

,  Common 

.  ?B  Reby  Gnd 

,  Common 

.Disc .R-02 

.Disc . R-10 

.Disc . R-ll 

.Disc .R-12 

.DR-00 

.DR-01 

.DR-02 

.DR-03 

,DR-0U 

.DR-05 

.DR-06 

. DR-0  7 

.DR-08 

.DR-09 

.DR-10 

.DR-11 

.DR-12 

.DR-13 

.DR-lk 

.  DR-15 

.DR-lo 

.DR-17 

.  NC 

.  DR-18 

.DR-19 

.  DR-20 

.DR-21 

.DR-22 

.DR-23 


TB3- 
TB*+- 


U0 

•Ul 

1*3 

1+1+ 


NC 
TBl+- 
ont . 1 . . 

NC 
ont.2.  TBU-1+7 

NC 

e  Cont  .  . 

.  .   TB7- 

Cont.  TB 

T37- 

TBU- 

NC 

TB5 


•1*5 

1*6 


TB1+-1+8 

79 

1*-1*9 

80 

50 


T36 


■51 
52 
53 
51* 
55 
56 
57 
58 

59 
60 
•6i 
62 
63 
61* 

65 
66 
61 
63 
69 
70 

TB7-71 
72 
NC 

TB7-73 
lh 
75 
76 
77 
78 


ADC 

Wiring 

Assignments 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

Ground 

NC 

Ground 

NC 

NC 

DSR-21 

In  o 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

DSR 
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Athena 

Athena 

J80300 

Plug-Pin 

Logic              Terminal    Board 

Number 

Name                      Assignments 

P10105-1 

. .BDR-00 

TB0-1 

.  .BDR-01 

2 

3 

.  . .BDR-02 

It 

.  .BDR-03 

4 

5 

.  .BDR-04 

5 

6 

..BDR-05 

6 

7 

..BDR-06 

7 

S 

. .BDR-07 

8 

9 

.  .3DR-08 

9 

10 

..BDR-09 

10 

11 

.  .BDR-10 

.       TB1-1 

12 

. .BDR-11 

2 

13 

.  .BDR-12 

3 

14 

..BDR-13 

4 

15 

.  .BDR-14 

5 

16 

.  .BDR-15 

6 

IT 

.  .BDR-lo 

7 

18 

..BDR-17 

8 

19 

.  .3DR-18 

9 

20 

.  .BDR-19       •  • 

10 

21 

..BDR-20 

.       T32-21 

22 

..BBR-21       .. 

22 

23 

.  .None            . . 

NC 

2U 

.  .BDR-22 

.       TB2-23       • 

25 

.  .Pitch    Data    1 .            2k       . 

26 

.  .Pitch    Data 

2.            25 

27 

.  .Pitch    Data 

3.            26 

28 

.  .Pitch    Data 

4.            27 

29 

.  .Pitch    Sign 

23       . 

30 

. . Yaw    Data    1 

29       • 

31 

. . Yaw    Data    2 

30       . 

32 

.  .Yaw    Data    3 

.  ..T33-31 

33 

. . Yaw    Data    k 

32       . 

3k 

. . Yaw    Sign 

33       . 

35 

. .Ace .R-00 

34         . 

36    . 

. .Ace .R-01 

35       • 

37    . 

. .Ace . R-02 

36       . 

38    . 

. .Ace . R-0  3 

37      . 

-     39    . 

. .Ace .R-0  4 

38       . 

ko    . 

. .Ace . R-0  5 

39       • 

ADC 

W  i  r  i  n  g 

Assignments 


P36-2  9 
-27 

25 

23 
21 
19 

17 

15 

13 

11 

9 

7 

5 

3 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

P35-13 

P35-11 

P35-9 

P35-7 

NC 

NC 

NC 

NC 

NC 

NC 

DSR-4 

DSR-5 
DSR-6 

DSR-7 
DSR-8 

DSR-9 


Fig. A. 9.1.  Wiring    assignments    for    J-Box    J80300    located 
under    floor    of   Athena    Computor 


(Continued  Fig. A. 9.) 
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Athena 
Plug-Pin 
N  umber 


P10205-30 
31 
32 
33 
31* 
35 
36 

37 

2  5 

39 

ko 

41 
42 

43 
44 
45 
U6 
47 
43 
49 
50 
51 
52 
53 
54 

55 
56 

57 
58 
59 
6o 
6i 
62 


Athena 

Logic 

N  ame 


DE-24 
DR-25 
DR-26 

DR-27 
DR-28 
Displ 
TI  Di 


TI  Di 
TI  Di 
NC 

TI  Di 


ay  Common 
splay-00 . 
-01. 

02. 

03. 

04. 

05. 
06. 

07. 
splay-08 . 
splay-09 • 

splay-10 . 
1 1 . 


J  80  300 
Terminal  Board 
Aasi  gnments 

TB8-81 
82 
83 
84 

55 

55 

87 

5  5 

89 

90 
.TB9-91 
92 
93 
94 
95 
96 

97 


•  TB9- 

•  TB9- 

.  .  NC 
.TB9- 


Burno 

Groun 

P3  Re 

Groun 

DISC 

Groun 


12. 

13. 
14. 
15. 
16. 

17. 

18. 

19. 
ut 
d 

lay  Groun 
d 

R-05 
d 


93 

99 
100 
TB10-101 
102 
103 
104 

105 


106 


.  .  NC 
.  .  NC 
d.TBl 
.  .  NC 
.TB10 
.  .  NC 


0-107 
-108 


ADC 

w 

iring 

Acs 

ignments 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

NC 

Fig. A. 9.1.    Wiring    assignments    for    J-Box    J30300    located 
under    floor    of   Athena    Computor 


76 


DSR-1 

Record 

2 

Step 

3 

IR  Gap 

1* 

IBM  Bit  1 

5 

IBM  Bit  2 

6 

2     it  it 

7 

IBM  Bit;  8 

8 

IBM  Bit  A 

9 

IBM  Bit  B 

BSR-35   . 

Ground 

P301-G 

H 
A 
B 

C 

D 

? 
M 


Note 


DSR-10  and  DSR-13  are  shorted 
together,  all  other  pins  on  both 
plug  are  not  connected. 


Fig. a. 9. 2.  Digital  Stepping  Recorder  ( DSR) -to-Athena  Cable 


All  pins  of  P300  A  and  B  are  wired  "one-to-one" 
to  ?35  and  P36 ,  respectively,  except  pins 
P300  B-ko    to  P300  B-5h    which  are  not  connected. 

Example:   P300B-20  P36-20 

P300A-37  P35-37 


Fig. A. 10.    Athena-to-ADC  Cable 
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C  o  mp  ut o r  Console  Register  Bit  Position  Function 

steering  register  bits  00-03         multiplexer  binary  address 
discrete  register  bit  02  tape  step-record  command 

discrete  register  bit  10  ei  d  of  record  cor.irr.and 

discrete  register  bit  11  address  set  command 

acceleration  register  bits  00-15     fourteen  bits  ADC  output 
Fig. A. 11   Listing  of  input  and  output  functions  and  their 
corresponding  Athena  computer  registers 

A. 5 -2   One  Channel  Sample  50  Tines  Per  Second  Program 

This  program  enables  the  sampling  of  one  analog  data 
channel  approximately  50  times  in  one  second  and  converting 
to  digital  form.   The  program  in  machine  language  with  some 
general  comments  is  given  in  Fig.  A.l^.   The  following 
sections  explain  the  subroutines  in  the  program.   Again, 
all  addresses,  core  locations,  instructions,  and  constants 
are  given  in  octal  code. 

A. 5. 3   Sample  Counting  Subroutine 

The  sample  counting  subroutine  controls  the  number 
of  samples  taken  and  indirectly  the  length  of  the  analog 
signal  that  can  be  sampled.   The  subroutine  extends  from 
address  02000-02006  (See  Fig.A.lU).   The  first  instruction 
02^007  at  address  02000  initializes  the  control  console 
and  is  really  not  part  of  the  counting  program.   Instruction 
06025o  at  address  02001  begins  the  counting  subroutine. 
It  clears  the  accumulator  by  transmitting  0000  0000  from 
Core  Location  (CL)  256.   Instruction  061+255  at  address 
02002  a-ds  0000  0001  stored  in  CL  255  to  the  cleared 
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Fig.  A. 12.   Connectors  for  wiring  harness  from  ADC  rack  to 
Athena  Computer  and  wiring  harness  from  the  DSR 

to  Athena 


Vj 


Fig.    A. 13.       J-Box    80jC0    located    under    floor    of   Athena 

Computer 
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ADDRESS      INSTRUCTION  COMMENTS 


02C00  024007 

02001  0c025o 

02002  064255  (counting  program) 

02003  00U256 

02004  066252 

02005  222007  sign  jump 

020  0b  20  0  331 jump  to  end  of  run  program 

02007  C20001 

02010  024012  sets  multiplexer  address 

02011  02^007 

02012  006031  stores  binary 

02013  022001 places  channel  number 

02014  024010 

02015  004100  1st  step-record 

02016  004100  generated 

02017  024007  


02020  202175  time  delay 

generated 


02175  062031 

02176  lOoOlo  complements  sample  and 

clears  accumulator 

02206  060027 

022C7  104004  shifts  1st  four  bits 

to  accumulator 

02212  004040 

02213  02204C places  bits  at  tape  input 

02214  024010 

02215  00410C  2nd  step-record 

02216  004100  pulse  generated 

02217  024007 

02220  202406  time  delay 

generated 

02406  104006 

6  more  bits  shifted  and 

02412  004040  placed  on  tape  input 

02413  022C40 

02414  02401C 


02415  004100  3rd  step-record 

02416  004100  pulse  generated 

02417  024007 


02420  202o06  time  delay 

gene  rate  d 

02606  lolooo 

02612  004040  last    6    bits    shifted    and 

02613  022040  placed    at    tape    input 

u2.<14  024010 

02615  004100  4th    step-recc  ■ 

02ol6  004100  nuice 


1 


.tea 


02ol?  024007 

02620  202000  timed  jump  to  beginning 

of  program 

Fig. A. 14.   50  Sample/sec  one  channel  sample  program 
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accumulator.   Instruction  0 0 U2 56  stores  the  accumulator 
in  CL  256.   Instruction  066252  at  address  02004  subtracts 
the  number  in  CL  252  from  the  number  in  the  accumulator. 
The  number  in  CL  252  is  equal  to  N+l  ,  where  N  is  the  number 
of  required  samples. 

Instruction  22200?  at  address  02005  is  a  sign  jump. 
If  the  quantity  in  the  accumulator  is  negative  and  not 
equal  to  zero,  the  next  instruction  executed  is  a.z    address 
02007.   This  is  the  case  that  exists  until  N  samples  have 
been  converted.   On  the  N+l  sample  the  number  in  the 
accumulator  is  equal  to  zero.   When  this  happens  the  sign 
jump  instruction  causes  the  next  instruction  in  sequence 
to  be  executed.   This  instruction  is  200331  at  address 
02006  which  is  an  unconditional  jump  to  the  instruction 
at  address  331.   At  address  331  the  End  of  Run  program 
is  stored. 

Assuming  that  a  negative  number  is  in  the  accumulator 
at  the  execution  of  the  sign  jump  instruction  (i.e.,  less 
than  N  samples  have  been  converted),  instruction  020001 
at  address  02007  is  executed.   This  instruction  transmits 
to  the  steering  register  the  sampled  channel  number  stored 
in  CL  1.   The  number  in  CL  1  is  7776  7777-   Only  the  low 
order  bits  are  transmitted  to  the  steering  register. 
True  bits  are  placed  on  steering  register  positions  00-05. 
The  pir.s  connected  to  these  positions  have  the  opposite 
sense  of  the  control  console  register  positions.   Therefore, 
no  voltage  appears  on  output  pins  corresponding  to  positions 
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00-05.   The  steering  register  is  vired  to  the  multiplexer 
binary  address  input  (Fig. A. 9),  thereby  setting  the  analog 
channel  to  be  sampled  to  zero. 

Instruction  C2UC12  at  address  0210  begins  the  required 
address  set  pulse.   This  pulse  must  be  at  least  HO  usee. 
long  and  equal  to  -oV  for  the  digital  stepping  recorder 
to  start  izs    step-record  process.   Instruction  02^012 
transmits  2000  stored  in  low  order  core  location  12  to  the 
descrete  register.   This  turns  off  position  11,  the  address 
set  command  bit.   This  bit  exhibits  a  reverse  sense  at  its 
output  pin  connected  to  the  multiplexer  (Fig. A. 4).   Therefore 
when  bit  11  is  turned  off  -6V  appears  at  the  multiplexer 
input.   Forty  usee. later,  the  next  instruction  02lt007  at 
address  02011  turns  on  bit  11  turning  off  the  -oV  at  the 
multiplexer  input. 

The  ADC  completes  digitization  39  usee. after  the  address 
set  pulse  is  initiated.   Therefore,  the  next  instruction 
executed  can  store  the  Ik   bit  ADC  binary  output  in  core. 
These  1 4  bits  are  shown  connected  to  Athena  in  Fig.A-.3. 
Instruction  006031  at  address  02012  stores  the  Ik    bits 
is  CL  31  - 

The  first  step  in  writing  the  magnetic  tape  is  to 
record  the  channel  number  on  the  first  tape  character  as 
shown  in  section  A. 6.   Instruction  022001  at  address  02013 
transfers  the  channel  number  to  the  acceleration  bit 
positions  00-05-   Again  the  corresponding  output  pins 
connected  to  the  Digital  tape  deck's  six  inputs  exhibit 
reverse  sense.   Thus,  the  complement  of  the  channel  number 
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is  transmitted  to  the  acceleration  register. 

A  step-record  pulse  of  -6V  and  longer  than  50  usee  . 
is  required  to  record  the  first  tape  character.   Instruction 
021+010  at  address  02ltll  transmits  000U  in  CL  10  to  the 
descrete  register.   Bit  position  02  of  the  discrete  register 
is  turned  on.   The  corresponding  output  pin  connected 
to  the  s^ep-record  input  is  of  the  correct  sense.   The 
length  of  this  pulse  is  increased  by  the  two  dummy  instructions 
00i+100  and  001+100  at  addresses  02015  and  02016  .   These 
instructions  store  the  accumulator's  content  in  CL  100 
which  in  no  way  affects  the  program.   Instruction  02^007 
s.z    address  02017  turns  of  discrete  register  pin  02  and  the 
-oV  applied  to  the  step-record  input.   The  first  tape 
character  has  now  been  written. 

The  next  instruction  is  an  unconditional  oumP  required 
for  the  synchronous  operation  of  the  digital  tape  deck's 
stepping  motor.   The  jump  instruction  is  202175  at  address 
02020.   The  next  instruction  to  be  executed  is  at  address 
02175  and  begins  the  process  of  writing  the  second  tape 
character.   Instruction  0o2031  a.z    address  02175  causes 
the  complement  of  the  1*+  bit  sample  store  in  CL  31  to  be 
transmitted  to  the  accumulator.   The  complement  is  trans- 
mitted to  account  for  the  acceleration  reverse  sense  pins 
connected  to  the  digital  tape  deck's  six  input  connectors. 

The  \U    sample  bits  are  shifted  left  into  the  quotient 
register  14  positions  so  that  the  writing  process  can  begin. 
Instruction  lOoClo  at  address  02176  does  the  shifting. 
This  shifting  process  takes  eight  instruction  times  so  that 
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the  next  instruction  is  at  address  022C6.   The  instruction 
is  O0OO27  and  transmits  0000  from  low  order  core  location 
27  to  the  accumulator  which  clears  the  accumulator  for  the 
sample  writing  process. 

The  14  complemented  sample  bits  are  shifted  four 
positions  to  the  left  by  instruction  104004  at  address 
02207-   This  peaces  the  first  four  of  the  negative  sample 
bits  on  accumulator  bit  positions  00,01,02,  and  03. 
Accumulator  bit  positions  04  and  05  are  of  no  consequence 
in  this  character  due  to  the  IBM  360/50  decoding  program 
blanking  them  out.   The  contents  of  the  accumulator  are 
stored  in  CL  40  by  instruction  004040  at  address  02212. 
Then  the  consents  of  CL  40  are  transmitted  to  the  accel- 
eration register  by  instruction  022040  at  address  02213. 
This  places  the  bits  of  the  second  tape  character  at  the 
proper  tape  deck  input  pins  and  of  correct  polarity. 

The  step-record  pulse  is  generated  by  instructions 
at  addresses  02214,  02215,  02216  ,  and  02217  as  previously 
explained. 

A  5ms  delay  for  synchonous  tape  drive  is  obtained 
by  the  jump  instruction  202406  at  address  02220.   In- 
struction 104006  at  address  02406  shifts  the  next  six 
complemented  sample  bits  from  the  quotient  register  into 
accumulator  bit  positions  00-05  •   Instruction  004040  at 
address  02412  stores  the  accumulator  conten-s  in  CL  40. 
Instruction  C22043  at  address  02413  transmits  the  contents 
of  CL  40  to  the  acceleration  register  placing  the  next 
six  complemented  sample  bits  on  acceleration 
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bit  positions  C0-05-  This  places  the  proper  tape  inputs 
for  generating  the  third  tape  character  at  the  tape  deck 
input.  A  step-record  pulse  is  generated  as  described  by 
instructions  at  addresses  02^14,  02^15,  02*116  ,  and  02*117. 

Instruction  202o06  at  address  02^20  is  a  jump  required 
for  the  synchronous  drive  of  the  tape  deck.   At  address 
026o6,  instruction  10ir006  shifts  the  remaining  complemented 
sample  bits  from  the  quotient  register  into  the  accumulator. 
These  four  bits  are  in  accumulator  bit  positions  CO,  01, 
02,  and  03.   Bits  in  pos.tions  Ck    and  05  are  ignored  by 
the  IBM  3oO/50  decoding  program.   The  accumulator  contents 
are  stored  in  CL  kO    by  instruction  001+040  at  address 
02ol2.   From  CL  10  the  contents  are  transmitted  to  the 
acceleration  register  such  that  the  proper  complemented 
sample  bits  are  placed  at  the  digital  recorder  input 
connections  for  writing  the  fourth  and  last  character. 
The  step-record  pulse  is  generated  by  instructions  at 
&26ll+,  02615,  026l6,  and  02617. 

Instruction  202000  at  address  02o20  is  an  unconditional 
jump  to  the  beginning  of  the  sample  program  and  starts 
the  process  over  again. 

A.  5.^-   End  of  Run  Subroutine 

The  end  of  run  subroutine  places  a  number  corres- 
ponding to  the  trial  number  placed  on  the  Athena  computer 
constant  register.   In  addition,  the  next  three  characters 
are  also  -written  with  the  run  number.   This  allows  the 
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360/50  computer's  decoding  program  to  treat  the  trail 
number  exactly  as  if  it  was  a  channel  number.   The  remain- 
ing characters  are  treated  as  if  they  were  a  sample  and 
must  be  ignored  by  the  programmer. 

The  run  number  program  was  written  on  drum  number 
three  in  such  a  fashion  so  that  the  tape  drive  pulses  are 
5ms  apart.   In  addition,  an  interrecord  gap  command  is 
generated  5ms  after  the  run  number  has  been  written  and 
sent  to  the  digital  tape  deck.   A  listing  of  this  program 
is  give.,  in  Fig.  A. 15  along  with  some  comments. 

A. 5. 5   Load  Constant  Subroutine 

The  load  constant  routine  must  be  run  before  any 
attempt  at  data  conversion  is  begun.   This  program  ini- 
tializes all  core  locations  used  in  the  counting  routine 
and  control  bits.   It  is  located  on  drum  number  zero 
between  addresses  00li-l  and  00211.   The  load  constant 
program  is  listed  in  Fig.  A.l6. 

A. 6   Definition  of  Binary  Data  Magnetic  Tape  Bit 
Structure 

The  seven  track  2^00  tape  units  at  the  Kansas  State 
University  computing  center  have  an  optional  feature  which 
allows  the  writing  and  reading  cf  binary  data  other  than 
Binary  Coded  Decimal  (BCD).   This  is  the  feature  -chat 
allows  the  ADC  digital  output  to  be  recorded  on  magnetic 
tape  with  the  DS3  _23C  stepping  recorder  and  ; ubseauently 
entered  into  the  IBM  3oO/50  core  for  translation  to  a  usable 
form. 
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ADDRESS 


UCTION 


COMMENTS 


06lTo 
061 77 
06200 
06201 
06202 
06203 
0620U 
06205 
06206 
06207 
06210 
06211 
06212 
06U06 
061*07 
061(10 
06H11 
061(12 

06606 
06607 
06610 
06611 
06612 

07006 
07007 
07010 
07011 
07012 

07206 
07207 
07210 
07211 


060027' 
001+256 

060027" 

CO 20 70 

001+001+" 

06200U 

0  0  h  0  0  H 

02200**" 

021+010" 

001*100 

001*100 

021*007 

2061*06" 

021*010" 

001*100 
001*100 

021+007 
2066O0" 

02l*C10" 

001*100 

OOllOO 

021*007 

207006" 

02U010" 

001*100 

001*100 

021*007 
207206 

01.7005 
060027 
017005 

32001*0 


:.  r 

TO 


Ltialize: 

t_ne  to 


counting 
zero 


cle  ars 
enters 


accumulator 
run  number 


complements  run  number 


places  run  number  at  tape  input 


1st  step -re cord 
pulse  generated 


timing    jump    delay 


2nd  step-record 
pulse  generated 


timing  j  umb 
delay 


3rd  step-record 

pulse  generated 


1 1  mi  n  g 
delay 


urn 


1+th    step-record 
pulse    generated 


timing 

jump 

delay 

interre  cord 

ga 

.p 

pulse 

generate 

c 

manual    stoo 


A. 15.   Enter  run  number  program 
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ADDRESS 

INSTRUCT 

l4l 

127776 

142 

147777 

11*3 

004001 

Ikk 

127775 

1U5 

147776 

146 

004002 

1U7 

12  7' 

150 

147775 

151 

004003 

152 

126000 

15  3 

140000 

154 

004007 

155 

126004 

156 

140000 

157 

oo4oio 

160 

124000 

161 

140000 

162 

004C11 

163 

122000 

164 

1-;0C00 

165 

004012 

166 

120000 

167 

140000 

170 

004027 

171 

121751 

172 

140000 

173 

004254 

nh 

120001 

175 

140000 

176 

004255 

177 

12000C 

2C0 

140000 

201 

004256 

202 

120004 

203 

140000 

204 

004004 

205 

125671 

206 

140000 

207 

004252 

210 

320004 

211 

200040 

Fig. A. 16 

COMMENTS 


loads  CL  1 


s  CL  2 


loads  CL  3 


loads  CL  7 


loads  CL  10 


loads  CL  11 


loads  CL  12 


loads  CL  27 


1751o  equals  N+l 

samples  (controls  length 
of  3  channel  program  can  be 
varied  as  required) 


loads  CL  25? 


loads  CL  2?b 


loads  CL  4 


56  71o  equals  N+l 

(can  be  changed  to  control 
length  of  5  0  s / s program) 


manual  stop 
jump  to  00040 


Program  to  load  constants 
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The  binary  conversion  feature  is  set  on  by  a  job 
control  card.   This  forces  the  seven  track  BCD  translator 
off  and  the  tape  parity  check-,  to  odd.   Therefore,  the  seven 
track  tape  is  written  with  odd  parity. 

The  binary  code  chosen  takes  advantage  of  the  fact 
that  four  tape  characters  are  placed  into  three  bytes  of 
the  computor  core  (IBM,  Ref.ll).   The  code  was  designed 
so  it  included  the  following  data, 

1.  multiplexer  address  of  channel  being  sampled 

2.  ill  bit  sarr.pled  data 

3.  run  number. 

Fig. A. IT  illustrates  how  two  samples  and  their  corresponding 
channel  numbers  appear  on  the  digital  magnetic  tape.   The 
run  number  is  written  in  the  address  location  after  the 
required  number  of  samples  of  data  have  been  converted. 
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Fig. A. 17      Digital    tape    bit    structure    definition 
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A. 6.1   IBM  3o0/50  Assembler  Language  Decoding  Program 
cc  n  track  tape  unit  converter  places  the  tape 
characters  into  the  computer  core  as  shown  in  Fig.A.lS 
(parity  bits  are  ignored). 


XXA4A3A2A1XX  D1D2D3DltD5D6DTD8  ^D^D^D^D^D^XX 
012  3^56 7  0123456  7   01   2   3  U      5   67 

Core  Locations 

Fig. A. 18.   Address  and  sample  bits  as  they  appear  in  core 


Tho  data  is  in  a  modified  negative  radix  binary 
notation  with  a  1  bit  positive  and  a  0  bit  negative. 
The . 360/50  assembler  language  program  was  written  to 
separate  the  data  from  unused  bits  and  write  the  uata  en 
tape  in  Fortran  acceptable  format  10(13,15). 

The  data  translation  from  negative  radix  is  accomp- 
lished by  manipulation  of  the  equation  for  measured  voltage 
reconstruction  as  follows 


Via  =  E 


1  3 


f 


B 

I    Bn) 


-I 


n-o   2n 


Vm   =    -f|if  P         Bn2132-4      -213 


Vm    =    -~ 


fs 

13 


'    n  =  o 

L 


"nBn)       -213  I 
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The  summation  can  be  evaluated  in  base  10  by  treating 
the  data  as  if  it  were  a  positive  radix  (base)  2  number 
and  converting  to  base  10  ir.  the  normal  manner,  as  shown 
in  Fig. A. 19- 

Original  Bit   SQ  B1  B2  B3  B^  B$  Bg  000  B13 

13-n       13  12  11  10  9   8   7   ooo  0 

Modified  Bit  Value  213  2~2  211  210  29  2   27  ooc  2° 

Fig.  A. 19-   Conversion  of  negative  radix  2  number  to 

base  10 

In  order  to  avoid  roundoff  errors,  only  the  expression 
in  brackets  is  evaluated  and  is  written  as  a  fixed-point 
number  in  80  column  card  images  in  blocked  format  on  a 
magnetic  tape.   The  remainder  of  the  conversion  is  done  with 
the  user's  Fortran  program  which  must  convert  the  five 
digit  integer  sample  to  a  real  number  and  multiply  by  the 
scale  factor  10/1892. 

The  assembler  program  is  completely  general  and  can 
be  used  for  any  data.   Input  and  output  blocksize  and  tape 
densities  can  be  varied  using  the  360/50  job  control 
language  recompling  the  assembler  language  program.   See 
Fig.  A. 20  and  Fig.  A. 21  for  a  flow  chart  and  listing 
of  the  assembler  convertor  program,  respectively. 

A. 6.2   Notes  to  the  Programmer  Operator 

The  length  of  the  analog  signal  to  be  converted 
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Fig. A. 20.   Data  Conversion  Flowchart 


//ASMBL   EXEC  ASMFCLG.PARM.  ASM=  'LOAD* DECK  »l_  I  NECNT=5&' 

//ASM. SYS  IN   DD   * 

TIT'  E  l  ATHENA  DATA  CONVERTER  J.P.L." 

CrNVRT    START  X'E820«  THIS  ST<  .    MT  IS  OPTIONAL. 

ORG       SAVE  (  14*12) «•*  SAVE  THE  F      'ERS. 

11*15  LOAD  THE  BASE  REGISTER. 

USING  ORG, II  ,  ;  ABLISH  ADDRESSABILITY. 

ST  13*SAVE+4  RE  HIS  SAV  AREA  ADDRESS. 

LA  5*SAVE  GET  CUP   A  .  ,REA  ADDRESS. 

ST  5,8(0,13)  STORE  IT. 

LA  MAKE  IT  AVAILABLE  TO  CALLED  PROGRAMS. 

OPEN  ( INDCB»( INPUT ),OUTDC3* (OUTPUT) )  OPEN  DATA  SETS  FOR  PROC. 

PllTDATA   PUT  OUTDCE  EMPTY  A  BUFFER  AND  GET  A  NEW  ONE. 

LR  4.1  COPY  ADDRESS  OF  BUFFER. 

LA  5.10(0,0)  LOAD  10  INTO  REGISTER  5. 

IV I  0(1),C*  '  SET  UP  TO  

MVC  1  {  131  ♦  1 )  »0("l)  CLEAR  THE  NEW  BUFFER  TO  BLANKS" 

GFTDATA   G;:T  INDCB,IN-:-l  GET  AN  INPUT  RECORD. 

NC  :N.=X'rC03CFFFC,     ZERO  THE  INPUT  AREA  EXCEPT  GOOD  BITS. 


XR     2»2  CLEAR  REGISTER  2» 

2?  IN  LOAD  THE  DATA  INTO  IT. 

SRDL   2,18  SHIFT  CHANNEL  ADDR.  TO  RIGHT  END. 


SRL    3,16  SHIFT  DATA  TO  RIGHT  END. 

O/D    2*DW0RD  CONVERT  CHANNEL  TO  PACKED  DECIMAL. 

CI     DWCRD+7»X*  OF' FIX  THE  SIGN  FOR  UNPACK. 

US  I NG  CU f»  4  AD  D  RE  S  S  AB I L  IT  Y  FOR  T  HE  0  U  T  PU  T  RECORD  . 

UNPK   OUTCHANtDWORD       UNPACK  THE  CHANNEL  NUMBER. 

CVD    3, DWORD  CONVERT  DATA  TO  PACKED  DECIMAL. 


SP     DWOED,=P' 8192'       COMPENSATE  FOR  DATA  SIGN. 

UNPK   OUTDATA, DWORD       UNPACK  THE  DATA  INTO  OUTPUT  AREA. 

TM     BITS»X"lO>  CHECK  THE  SIGN  OF  THE  DATA. 


B2     FIXSIGN  BRANCH  IF  IT  IS  PLUS. 

MVI.    OUTSIGN»C'-'        OTHERWISE*  MOVE  IN  A  MINUS  SIGN. 
FIXSIGN   01     BITS»X*FO'  WIPE  OUT  THE  ORIGINAL  SIGN. 


A      4»  =  F«3'  INDEX  DOWN  THE  OUTPUT  AREA". 

JCT    5*GETDATA  LOOP  TEN  TIMES. 

5      PUTDATA  WRITE  OUT  THE  COMPLETED  RECORD. 


DROP   4 

quit    close  (outdcb*, indcs)    close  data  sets  and  write  last  rec. 
_l    13»save+4        reload  his  savearea  address. 


RETURN  (14,12!  RESTORE  THE  REGISTERS  AND  QUIT. 

INDCB     DC8    DSORG=PS»MACRF=GM»DDNAME=INPUT,RECFM=FB,EODAD=QUIT» 
DEVD^TA 


OUTDCB    DCB    DSORG  =  PS  *MACRF=PL  ,  DDNAME=CUTPUT  *RECFM=FB  » 

DEVD=TA 

_SAVE DC 1  8F'0' 

I N        DS 

DWORD     DS     ID 
LTORG 


OUT       DSECT 

DS  CL1 

OJJTCHAN   DS CL2 

OUTSIGN   DS  C'Ll 

3  T  )  VA   DS  CL4 

ORG_  OU 

"BITS      DS~~  CL1 

END  ORG 
/* 


Fig. A. 21.  Computer  Listing  of  Convertor  Program 
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es  the  time  the  sample  program  must  run.   The 
length  of  the  50  sample/sec  one  channel  program  is  controlled 
by  the  octal  number  stored  in  core  location  252.   The 

octal  number  occupies  both  the  high  and  lev  order  core 
positions.   In  Fig. A. 14,  5t>71->  equals  ^001.  converted 
samples.   This  is  stored  in  the  low  order  positions  by 
instruction  125671  at  address  00205.   The  high  order 
zeros  were  made  zero  ay    instruction  ll+OOOO.   This  program 
must  be  run  after  the  constants  are  entered  in  their  proper 
place  in  this  program  for  them  to  be  placed  in  core. 

Calculate  the  number  to  be  placed  in  CL  252  by  con- 
verting the  analog  signal  time  duration  to  samples  by 
the  conversion  factor  50  samples  =  1  sec.   Usually,  the 
number  of  samples  chosen  is  an  integral  multiple  of  10 
to  ease  readout  of  the  data  at  the  IBM  360/50.   Do  not 
forget  to  convert  the  decimal  number  to  octal  before 
entering  it  into  its  proper  place. 

It  is  a  good  idea  to  run  the  sample  program  once  with 
the  digital  tape  recorder  0?F  to  check  the  length  of  run 
and  to  initialize  the  computor. 

The  run  number  is  entered  on  the  first  two  registers 
of  the  constant  program  in  OCTAL  form. 

After  all  data  has  been  converted  punch  END  OF  FILE 
(EOF)  on  the  front  panel  of  the  digital  tape  recorder. 
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.'.NDIX  B 


It  is  desired  to  show  that  tl    .  urier  transform  of  Eq. 


5.24, 


°r..  (t)  =  h2(t)0{h1(t)*[h1(t)0<I.ff(T)]} 


(B.l) 


is  equal  to  , 


$rs(jw)  =  H2(-Jw)|Hi(ju)|'6*   (ju)   , 


(B.2) 


where  the  symbol  3  means  correlation  in  the  aperiodic  sense  and 
the  symbol  *  denotes  convolution. 


The  defining  equation  for  the  Fourier  transofr: 
f(t)  is, 


of  a  function 


F(ju)  =  /   f(t)e  jWtdt 


(3.3a) 


which  leads  to  the  equation, 


F(-jw)  =  / 


f(t)ejUitdt 


(B.3b) 


To  begin  let  g..  (t)  equal, 


h1<t)*[h1(t)e*ff(T)3  , 


(3.4) 


in  Eq .  3.1.   Eq.  B.l  written  with  the  substitution  indicated  by 
Eq .  3.4  is , 


.-;,  -  h2(t)0gl(t)   . 


(3.5) 


The    F    uric       transform    c  (-,)    is 
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OO  00 


(joi)    =    /      /      h     (v)g     (v+x)dve 

r  s  z  ■  _ 

—  OO        — OO 


v    .         -i  OJT   - 


/      h2(v)e-1       dvf         g 

—  OO  _0O 


(v+T)e"-j(V'T)dT       , 


which  by  comparison  with  Sqs.  B.3a  and  B.3b  means, 


*  0  (jw)  =  H,(-jw)G,  (jw) 


(B.6) 


The  Fourier  transform  G..  (jto)  becomes,  using  the  operations 
indicated  by  Eq.  B.4, 


G-(ju)  =  /  /  h  (a)g 2(t-a)dae"ja3tdt   , 


—  CO  — CO 


where  , 


82(t) 


=  h  (t)  0  $-,(T)       and, 


OO  CO 

G1(ju)  =  /  h1(a)e_-1a)ada/  g2(t-a)e_ja)(t-a)dt   , 


which    again    by    comparison    with    Eqs.     B.3    becomes, 


G,  oa)    =    H    (  jw)G,  (jw) 


(B.7) 


Combining  Eq.  3.7  with  Eq .  B.6  yields, 


^rs(jw)  =  H2(-ja))H1(jw)G2(jw)   . 


(B.8) 


57 


ACKNOWLEDGEMENT 


I  wish  to  express  my  appreciation  to  my  advisor,  Dr. 
Donald  Lenhert,  Drs.  W.  W.  Koepsel  and  M.  R.  Fedde  for  their 
advice,  understanding  and  material  contributions  during  the 
course  of  this  study. 

To  my  wife,  who  has  .served  as  an  able  assistant  and 
exhibited  confidence  in  me  throughout  this  work,  this  thesis 
is  dedicated. 


VITA 
CHARLES  H  .  II  CG  ITOWER 
Canii d;  tc  fo c    the  Degree  of 

MASTER  OF  SCIENCE 


Thesis:   The  Analysis  of  Nerve  and  Muscle  Electrical 
Waveforms  with  the  Fast  Fourier  Trans  form 
and  Correlation  Techniques 

Major  Field:   Electrical  Engineering 

Biographical: 

Personal  Data:   Born  in  Philadelphia, 

Pennsylvania,  July  1,  1944,  the  son  of 
Russell  L.  and  Margaret  P.  Hightower. 

Education:   Graduated  from  Canoga  Park 
High  School  in  1962;  received  the 
Bachelor  of  Science  Degree  from  Kansas 
State  University,  in  Electrical  Engineering, 
in  June,  1964;  completed  requirements  for 
the  Master  of  Science  Degree  in  June,  1968. 

Scholastic  Societies,  Professional  Organizations , 
ana  Honors:   Graduated  Cum  Laude  from  Kansas 
State  University  in  June,  1966;  elected  to 
membership  in  Phi  Kappa  Phi,  Eta  Kappa  Nu , 
and  Sigma  Tau ;  participated  in  Engineering 
Honors  Program  during  undergraduate  work; 
member  Institute  of  Electronic  and  Electrical 
Engineers  ;  completed  Engineer-in-Train:lng 
examination,  state  of  Kansas,  1966. 

.'rofessicn:.!  Experience:   Associate  Engineer  at 
Douglas  Aircraft  Company  working  en  Long  Tank 
Thor  guidance  system  from  June  to  September, 

associated 


THE  ANALYSIS  OF  NERVE  AND  MUSCLE  ELECTRICAL  WAVEFORMS  WITH 
.  E  FAST  FOURIER  TRANSFORM  AND  CORRELATION  TECHNIQUES 


by 


CHARLES  HOWARD  HIGHTOWER 
B.  S.,  Kansas  State  University,  1966 


AN  A3STRACT  OF  A  MASTER'S  THESIS 


submitted  in  partial  fulfillment  of  the 


requireaients  for  the  degree 


MASTER  OF  SCIENCE 


Department  of  Electrical  Engineering 


KANSAS  STATE  UNIVERSITY 
Manhattan,  Kansas 

1968 


ABSTRACT 


The  respiratory  system  common  to  a  large  number  of  living 
organisms  is  dependent  on  a  group  of  receptors  sensitive  to 
blood  carbon  dioxide  (CO-)  content  for  the  manipulation  of 
ventilation.   Collectively,  these  receptors  are  called  C0- 
chemoreceptors .   Their  specific  function  is  to  communicate  with 
the  central  nervous  system  and  notify  it  of  the  blood  C0- 
content.   The  central  nervous  system  through  a  complex  array 
of  nerve  bundles  then  stimulates  respiratory  muscles  to  perform 
the  operation  of  ventilation. 

Observation  of  the  electrical  activity  in  the  respiratory 
nerve  and  muscle  systems  yields  no  apparent  deterministic 
relationship  between  the  state  of  the  CO-  chemoreceptors  and 
the  electrical  activity  (which  is  mainly  stochastic  in  nature) . 

The  effect  of  the  CO-  chemoreceptors  on  respiratory 
electrical  activity  in  the  Gailus  Domesticus  was  investigated 
in  this  thesis  using  techniques  of  modern  communications 
theory  in  conjunction  with  data  processing  by  a  large  digital 
computer.   Results  of  this  analysis  indicated  that  the  effect 
of  the  CO-  chemoreceptors  on  respiratory  electrical  activity 
does  manifest  itself  in  two  deterministic  ways.   The  first 
descriptive  quantity  was  found  with  the  Fast  Fourier  Transform 
of  the  nerve  and  muscle  waveforms  to  be  the  bandwidth  of  the 
e!  -ctrical  waveform.   The  second  quantity  was  four...  to  be 

i   relative  power  content  of  the  electrical  wa\  -1 ;.s  as 

indicated  by  their  cumulative  power  spectra. 


The  behavior  of  these  two  descriptors  suggested  that 
the  C09  chernoreceptor  relationship  with  the  respiratory 
system  is  best  characterized  by  two  distinct  mechanisms . 
One  mechanism  predominant  during  phasic  respiration  and  the 
other  during  apnea. 

An  idealized  mathematical  model  of  the  respiratory 
system  was  constructed  and  analyzed,  showing  that  the 
Fourier  transform  of  the  nerve-muscle  crosscorrelation 
function  is  an  indicator  of  the  accuracy  of  such  a  model. 

In  addition,  a  brief  examination  of  aliasing  duo  to 
waveform  digitization  at  a  rate  below  the  Nyquist  frequency 
was  made.   From  the  results  it  was  recommended  that  the 
waveform  digitization  rate  be  increased  by  a  factor  of  from 
two  to  four  to  lessen  possible  aliasing  perturbations. 
Other  recommendations  were  made  concerning  the  experimental 
data  recording  equipment.   Also,  further  areas  of  research 
are  discussed  based  on  the  findings  of  this  thesis. 


